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1 . INTRODUCTION 


1 . 1 Program  Obj  ective 

The  ultimate  goal  of  the  work  reported  herein  is  to  develop  an  ob- 
jective method  for  the  short  range  forecast  of  storm  development  and 
motion.  The  initial  step  in  this  program  is  to  devise  a set  of  para- 
meters for  the  characterization  of  weather  radar  data  to  efficiently 
represent  the  essential  information  obtained  by  a radar  without  requiring 
extensive  storage  capacity  to  handle  unprocessed  data.  In  this  report 
we  consider  techniques  for  the  representation  of  the  reflectivity  and 
Doppler  velocity  fields  generated  by  a single  weather  radar.  The  reflec- 
tivity data  are  considered  both  alone  and  in  conjunction  with  simultaneously 
obtained  Doppler  data. 

A computer  program  was  developed  to  process  Doppler  weather  radar 
data  to  obtain  the  required  parameters.  The  program  detects  small  con- 
vective cells  and  larger  echo  regions  and  computes  a series  of  attributes 
for  each.  The  program  represents  the  first  step  in  the  development  of  an 
objective  procedure  for  the  automatic  processing  of  weather  radar  data 
for  use  in  the  short  range  forecast  of  storm  development  and  motion. 

1 . 2 Summary 

The  recommended  parameterization  of  radar  data  is  based  upon  the 
use  of  small  convective  cells  to  represent  the  basic  architecture  of  a 
storm  system.  Convective  cells  are  readily  apparent  in  isolated  showers, 
clusters  of  showers,  and  squall  lines.  They  are  also  evident  as  im- 
bedded structures  in  the  rain  bands  associated  with  widespread  rain. 

Crane  (1976)  found  that  small  convective  cells  were  stable  entities 
which  could  be  reliably  identified  on  successive  scans  and  tracked  from 
scan  to  scan. 

The  small  cells  are  characterized  by  a set  of  attributes;  intensity, 
area,  height,  age,  stage  of  development,  associated  low  level  convergence 
(radial  shear) , associated  vorticity  (tangential  shear) , and  propagation 
velocity.  Additional  parameters  are  obtained  to  characterize  the  cells 
within  a larger  precipitation  (echo)  region.  These  pareimeters  include 
cell  spacings  and  relative  orientation,  number  of  cells  within  a precip- 
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itation  region,  relative  cell  motion,  and  the  motion  of  the  cells  rela- 
tive to  the  motion  of  the  centroid  of  the  encompassing  echo  region.  The 
mean  radial  velocity  data  are  also  processed  to  estimate  the  mean  wind  pro- 
file (environmental)  and  to  identify  local  maxima  in  tangential  shear. 

The  shear  maxima  may  not  be  coincident  with  a single  cell  but  may  occur 
within  a ''ell  cluster.  The  location  of  the  shear  maxima  relative  to  the 
location  of  neighboring  cells  is  also  used  to  characterize  the  Doppler 
velocity  field. 

A previous  analysis  of  available  aircraft  observations  of  velocity 
fluctuations  and  of  radar  observations  of  tangential  shear  and  Doppler 
velocity  variance  by  Crane  (1976)  had  shown  that  the  velocity  variance 
was  primarily  caused  by  shear  within  the  radar  sampling  volume.  The  var- 
iance data  therefore  may  not  be  useful  for  the  estimation  of  the  intensity 
of  turbulence  as  described  by  an  eddy  dissipation  rate.  For  this  reason, 
the  mean  Doppler  velocity  estimates  provide  the  principal  data  to  be  pro- 
cessed. These  data  are  used  to  develop  radial  and  tangential  shear  esti- 
mates for  association  with  detected  cells  and  to  locate  tangential  shear 
maxima  not  associated  with  a cell.  Variance  data  are  used  to  mark  regions 
with  larger  than  normal  velocity  fluctuations  that  should  not  be  included 
in  estimates  of  the  environmental  wind.  Local  maxima  in  velocity  variance 
are  also  detected  for  comparison  with  the  attributes  of  local  maxima  in 
tangential  shear  to  test  the  hypothesis  that  the  major  contribution  to  the 
observed  variance  is  due  to  larger  scale  shear  rather  than  turbulence. 

The  data  parameterization  reduces  the  amount  of  data  required  to 
represent  the  initial  radar  observations.  Each  volume  scan  is  repre- 
sented by  detected  cells,  by  larger  echo  regions,  by  tangential  shear 
maxima,  and  by  the  attributes  of  the  cells  and  larger  echo  regions. 
Additional  information  is  provided  to  describe  the  spatial  organization 
of  the  cells.  These  data  will  be  used  in  the  forecast  of  cell  propa- 
gation as  defined  by  their  development  and  motion.  The  cells  and  their 
attributes  are  important  for  the  identification  of  severe  weather  and 
aircraft  hazards,  however,  they  do  not  represent  the  total  production 
of  precipitation  within  the  echo  envelope  surrounding  the  cells.  Addi- 
tional data  will  be  provided  to  represent  the  equivalent  precipitation 
depth  (accumulation)  within  a larger  echo  region. 


-5- 


1.3  Software  Development 


The  goal  of  this  contract  with  the  Air  Force  Geophysics  Laboratory 
(AFGL)  is  to  provide  computer  software  to  obtain  the  parameters  required 
to  represent  weather  radar  data.  The  radar  used  to  provide  the  data  is 
the  C-band  weather  radar  operated  by  the  Weather  Radar  Branch  of  AFGL  at 
Sudbury,  Massachusetts.  The  computer  programs  were  prepared  for  the 
CDC-6600  at  AFGL. 

Table  1 provides  a list  of  the  cell  and  echo  area  attributes  recom- 
mended as  important  for  the  efficient  representation  of  the  radar  data. 
Due  to  the  limited  duration  of  this  contract  it  was  not  possible  to  pro- 
vide software  to  obtain  all  the  attributes  on  the  list.  The  attributes 
identified  by  asterisks  are  calculated  by  the  first  generation  computer 
program  developed  under  this  contract.  These  attributes  describe  radar 
data  obtained  on  a single  azimuth  scan.  Algorithms  exist  to  combine  data 
from  a series  of  azimuth  scans  within  an  elevation  scan  (Crane,  1976)  but 
were  not  included  in  the  first  generation  computer  program. 

Although  cell  tracking  algorithms  are  also  available  (Crane,  1976) 
they  were  not  included  in  the  first  generation  program  package.  The 
development  and  fine  tuning  of  the  tracking  algorithms  require  experience 
with  the  cell  detection  program  under  a number  of  different  environmental 
conditions  such  as  isolated  showers,  squall  lines,  and  widespread  rain. 
Neither  the  data  nor  the  time  were  available  to  process  the  required  data. 

1.4  Organization  of  the  Report 

A review  of  radar  data  processing  is  given  in  Section  2.  Reflecti- 
vity-based parameters  are  discussed  in  Section  3.  The  use  of  Doppler 
data  is  considered  in  Section  4.  The  computer  algorithms,  a description 
of  the  software  package  and  sample  results  are  given  in  Section  5. 

Section  6 summarizes  the  results  obtained  to  date. 


TABLE  1 


2.  BACKGROUND 


Although  weather  radar  data  have  been  operationally  available  for 
: ^ many  years,  they  have  not  been  used  in  routine  objective  forecast  proce- 

dures. Weather  radar  data  were  initially  displayed  as  echo-filled  or 
echo-free  regions  on  a plan  position  indicator  (PPl)  display.  The  data 
displays  were  useful  in  locating  precipitation  regions  and  providing 
forecast  verification,  but  they  were  not  useful  for  measuring  storm 
’ intensity  or  displaying  the  structure  of  the  storm.  Next,  reflectivity 

data  were  depicted  using  fixed  level  contours.  These  contours  provided 
a graphic  display  of  storm  intensity  and  structure  but  unfortunately  only 
a limited  number  of  contours  could  be  displayed  and  interpreted.  Recently, 
the  use  of  color  displays  has  increased  the  number  of  contours  that  can 
^ be  displayed.  Operator  interpretation  is  still  difficult  and  the  data 

I require  additional  processing  before  they  are  available  for  quantitative 

I objective  analysis. 

f Digitized  radar  data  are  required  for  objective  analysis.  Most 

\ current  research  radars  obtain  and  store  the  radar  data  in  a digital 

form,  and  operational  systems  are  being  improved  to  provide  digital 
data.  Attempts  have  been  made  to  use  digitized  fixed  contour  level  data 
for  the  objective  forecast  of  storm  motion.  Recently,  Elvander  (1976) 
reported  on  the  performance  evaluation  of  three  different  techniques  to 
estimate  and  forecast  echo  motion.  The  first  (oldest)  technique  used  a 
linear  least  squares  tracking  procedure  to  follow  the  centroids  of  echos 
defined  using  a fixed  reflectivity  level  contour  (Barclay  and  Wilk, 

1970;  Wilk  and  Gray,  1970) . The  predicted  location  of  an  echo  region 
was  estimated  by  extrapolation  along  a least  square  curve  fit  to  the 
previously  observed  echo  locations.  This  procedure  can  not  handle  storm 
development,  growth,  or  decay  - only  storm  translation. 

The  other  two  tracking  techniques  used  echo  velocity  estimates 
based  upon  correlation  analysis.  Correlation  analyses  have  been  used 
for  years  to  study  echo  characteristics  and  their  changes  (see  for 
example  Kessler  and  Russo,  1963).  Recently,  two  separate  correlation 
procedures  were  tried  to  automatically  derive  storm  motions:  correla- 
tions involving  only  isolated  echo  regions  (Duda  and  Blackmer,  1972; 
Blackmer  et  al.  1973)  and  correlations  using  the  entire  PPI  display 
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(Austin  and  Bellon,  1974).  The  first  provides  independent  velocity  esti- 
mates for  each  echo  region;  the  latter  uses  a single  velocity  for  all 
the  depicted  echo  areas.  Crane  (1976)  <^ound  that  the  echo  areas  propa- 
gate to  encompass  the  growth  and  decay  of  small  enclosed  convective 
cells.  The  small  cells  have  regular  tracks  but  cells  within  a larger, 
isolated  echo  region  may  move  in  slightly  different  directions.  The 
motions  of  the  larger  echo  regions  were  erratic  as  they  merged,  separated, 
and  changed  to  encompass  the  developing  cells. 

Elvander  found  that  the  objective  procedures  that  forecast  the 
motion  of  echo  centroids  defined  using  the  lowest  level  reflectivity 
data  worked  best  when  based  upon  velocity  estimates  generated  using 
correlation  techniques.  He  reported  that  the  least  squares  curve  fit 
approach  worked  best  when  the  echo  centroids  were  defined  using  vertically 
integrated  liquid  water  content  (VIL)  data.  Since  the  VIL  values  are 
largest  within  the  small  active  regions  of  convection,  the  VIL  results 
should  be  similar  to  those  reported  by  Crane  (1976)  when  the  echo  regions 
are  dominated  by  a single  intense  cell.  The  National  Hurricane  and 
Experimental  Meteorology  Laboratory  has  also  been  experimenting  with 
objective  echo  identification  and  tracking  procedures  (Ostlund,  1974; 
Wiggert  et  al,  1976).  They  initially  used  the  echo  centroid  tracking 
procedure  but  have  recently  abandoned  that  technique  to  use  a procedure 
that  tracks  reflectivity  maxima  or  peaks.  The  locations  of  the  peaks 
(reflectivity  maxima)  are  found  by  best  fitting  (correlating)  the  observed 
reflectivity  values  with  a number  of  two-dimensional  Gaussian  distri- 
butions. The  best  fit  Gaussian  distribution  are  used  to  identify  the 
peaks  within  an  echo  to  be  tracked.  This  procedure  was  devised  to 
improve  the  operation  of  their  program  when  splits  or  merges  occur. 

An  alternative  development  in  the  representation  of  reflectivity 
maxima  or  peaks  within  a larger  echo  region  is  the  use  of  small  cells 
defined  by  contours  a fixed  level  below  local  reflectivity  maxima  within 
larger  echo  regions  (Crane,  1976).  These  cells  are  defined  by  small 
reflectivity  changes  and  correspond  to  volumes  that  encompass  updraft 
regions  during  the  growth  stage  of  cell  development  and  encompass  down- 
draft  regions  during  the  mature  stage.  They  are  defined  on  a single  scan 
by  local  reflectivity  maxima  only  a few  dB  above  their  surroundings. 

The  local  concentrations  of  liquid  water  are  reliably  detected  throughout 
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the  active  stages  of  cell  development.  Single  identifiable  regions  of 
locally  increased  liquid  water  content  persist  from  scan  to  scan  for  dur- 
ations of  5 to  50  minutes.  The  small  cells  are  continuous  in  height  and 
display  smooth  regular  horizontal  motion. 

Doppler  velocity  observations  show  that  the  small  active  cells  are 
important  elements  in  organizing  deviations  in  the  flow  field  from  that 
of  the  surrounding  or  environmental  flow  pattern.  Reported  Doppler 
velocity  measurements  show  little  deviation  from  the  environmental  or 
background  winds  over  much  of  the  volume  enclosed  within  an  echo  region. 
Doppler  velocity  observations  near  the  small  cells  reveal  the  convergence 
patterns  required  to  feed  the  updrafts  and  respond  to  downdrafts . The 
data  also  reveal  mesoscale  cyclones  (and  anticyclones)  associated  with 
the  updraft  regions  and  with  secondary  flows  caused  by  a number  of 
closely  spaced  cells. 

Currently,  the  analysis  of  single  Doppler  radar  data  is  based  upon 
comparison  with  simplified  kinematic  models  for  the  flow  fields  of 
importance  to  severe  weather:  supercells,  tornadoes  and  low  level  gust 
fronts  (Donaldson,  1970;  Browning  and  Foote,  1976;  Burgess,  1976;  Brown 
and  Lemon,  1976;  Zmic  et  al,  1976).  The  identification  of  regions  of 
severe  weather  is  made  by  comparing  the  Doppler  observations  with 
signatures  representative  of  each  of  the  models.  The  Doppler  data 
provide  a measure  of  the  severity  of  the  weather  associated  with 
features  of  the  reflectivity  field.  The  reflectivity  data  in  turn 
provide  the  means  to  forecast  the  motion  of  the  active  regions  that  are 
probable  sites  of  severe  weather. 

Doppler  data  have  been  mainly  used  for  the  display  of  the  flow 
fields  within  an  echo  region  (especially  multiple  Doppler  radar  data)  and 
for  the  identification  of  severe  weather.  They  have  not  been  used  in  an 
objective  fashion  to  forecast  the  motion  of  severe  weather.  Initially, 
in  the  objective  analysis  algorithms  developed  under  this  contract,  the 
reflectivity  data  are  to  be  used  to  identify  cells  and,  using  observations 
on  successive  scans,  their  motion.  The  Doppler  data  will  be  used  for 
the  identification  of  regions  of  severe  weather  or  possible  hazard.  The 
data  will  be  processed  in  a manner  to  allow  ready  incorporation  of  addi- 
tional features  of  either  the  Doppler  or  reflectivity  fields  that  appear 
to  be  important  after  detailed  analysis  of  a large  set  of  data  using  the 
initial  processing  algorithms. 
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3.  OBJECTIVE  ANALYSIS  OF  REFLECTIVITY  DATA 


Objective  analysis  of  reflectivity  data  must  provide  information 
for  use  in  forecasting  the  location  and  development  of  severe  weather 
and  for  use  in  measuring  the  production  of  precipitation  and  the  resul- 
tant distribution  of  precipitation  on  the  ground.  The  analysis  algorithms 
developed  under  this  contract  include  two  types,  (1)  the  small  cell 
analysis  using  peak  reflectivity  reference  contour  levels,  and  (2)  larger 
echo  area  analysis  using  fixed  echo  contours.  The  former  is  recommended  i 

because  of  the  association  between  convectively  active  regions  and  severe 
weather  and  because  of  the  utility  of  the  small  cells  for  the  forecast 
of  pattern  development  and  motion.  The  latter  is  recommended  to  keep 
track  of  the  precipitation  produced  by  the  active  cells.  No  attempt  will 
be  made  to  partition  the  precipitation  within  an  echo  region  by  cell. 

3.1  Small  Cells 

The  use  of  objective  techniques  for  the  detection  of  small  convec- 
tive cells  was  developed  and  reported  by  Crane  (1976).  He  found  that  a 
small  cell  can  be  readily  detected  using  at  most  three  azimuth  scans; 
the  detection  probability  for  a single  scan  was  above  0.6  for  the  reflec- 
tivities greater  than  35  dBZ  and  for  three  scans  in  a volume  scan  sequence, 
the  probability  of  detection  increases  to  0.93.  The  detection  probability 
is  still  higher  for  a typical  volume  scan  with  a larger  number  of  azimuth 
scans  at  different  elevation  angles. 

The  small  cell  detection  procedure  developed  by  Crane  is  illustrated 
in  Figures  1 and  2.  Figure  1 shows  a hypothetical  echo  region  (lowest 
level  contour)  including  two  cells.  The  cells  are  identified  by  smaller 
contours  T*  units  below  their  enclosed  relative  maxima.  The  cell  areas 
are  the  shaded  regions  within  the  peak  referenced  contours.  Figure  2 
shows  the  application  of  this  detection  process  to  actual  radar  data  using 
a 2.5  dB  value  for  T.  The  outer  or  lowest  level  contours  on  this  figure 
have  a value  of  20  dBZ.  The  peak  values  are  above  50  dBZ.  The  data 
reveal  a wealth  of  detail  not  evident  using  a limited  number  of  fixed 
level  contours  separated  by  large  differences  in  reflectivity.  The  dis- 
play as  shown  in  Figure  2 is  quite  complex.  It  may  be  replaced  by  the 


display  in  Figure  3 with  little  loss  of  information.  On  this  figure  the 
cell  attribute,  peak  reflectivity,  is  listed  for  each  cell.  Other  attri- 
butes such  as  cell  area,  cell  height,  height  of  maximum  reflectivity, 
height  of  cell  base  for  first  echo,  or  any  other  measurable  associated 
with  the  reflectivity  of  Doppler  data  fields  may  be  calculated  and  dis- 
played or  recorded  for  each  of  the  cells. 

3.2  Larger  Echo  Areas 

The  small  cells  are  generally  contained  in  larger  multicell  echo 
regions.  Warner  (1976)  reported  that  all  the  hailstorms  he  observed  in 
Alberta,  Canada  occurred  as  small  cells  within  larger  echo  areas.  An 
analysis  of  his  reported  data  shows  that  single  isolated  cells  do  not 
develop  significantly  either  in  height  or  intensity.  Clusters  of  cells 
within  a single  envelope  defined  by  a low  level  reflectivity  contour 
(10-20  dBZ)  usually  exhibit  more  significant  development  growing  both 
higher  and  more  intense  than  the  single  isolated  cells  that  appeared  at 
the  same  time  on  the  same  day. 

The  apparent  cooperation  between  small  closely  spaced  cells  has 
been  reported  by  other  investigators.  Woodley  and  Simpson  (1972)  have 
reported  that  convective  cells  in  the  Florida  area  have  a higher  inten- 
sity and  produce  more  precipitation  after  they  merge  than  before.  They 
declare  mergers  when  echo  regions  defined  by  a fixed  25  dBZ  contour 
combine  to  form  a larger  multicell  echo  region.  Their  data  show  that 
the  environment  surrounding  each  cell  is  important.  The  number,  spacing, 
and  relative  orientation  of  the  cells  within  a single  echo  region  appear 
to  affect  the  development  of  the  small  cells.  These  data  must  be 
recorded  in  addition  to  the  attributes  of  each  cell.  They  are  attri- 
butes of  the  larger  echo  region. 

The  small  cells  are  considered  to  be  the  active  regions  of  convec- 
tion within  the  larger  echo  region.  The  larger  area  encompasses  precipi- 
tation resulting  from  the  transport  of  liquid  water  (and  ice  or  snow) 
and  water  vapor  out  from  the  updraft  regions.  The  transport  processes 
are  mainly  turbulent  - eddy  diffusion  or  advection  depending  upon  the 
scale  size  of  the  motion.  Microphysical  processes  continue  to  produce 
precipitation  within  the  larger  regions  about  each  cell  which  results  in 
precipitation  that  is  measurable  on  the  ground.  The  precipitation  is 
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apparently  carried  out  from  the  cells  by  the  environmental  or  background 
winds  as  it  settles  to  the  ground.  The  larger  echo  area  also  contains 
the  decaying  cells  that  remain  after  their  active  stages  have  been  com- 
pleted. The  total  precipitation  in  this  region  is  important  and  must  be 
obtained  from  the  radar  data.  The  transport  processes  are  complex  and 
it  does  not  appear  reasonable  to  attempt  to  identify  the  resultant  pre- 
cipitation with  particular  active  generating  regions. 

At  midlatitudes  the  larger  area  of  precipitation  surrounding  the 
active  cells  generally  consists  of  ice  and  snow  aloft  melting  to  form 
rain  below.  Care  must  be  taken  in  processing  the  data  to  exclude  measure- 
ments made  within  the  melting  region  or  bright  band.  Data  taken  at  the 
lowest  elevation  angle  will  be  processed  once  per  volume  scan  to  provide 
an  estimate  of  the  rain  rate  integrated  over  the  area  of  the  larger  echo 
region.  Rain  rate  estimates  made  on  a series  of  scans  (lowest  elevation 
from  each  volume  scan)  will  be  combined  to  estimate  the  averaged  accumu- 
lation of  precipitation  at  the  surface. 
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4.  OBJECTIVE  ANALYSIS  OF  DOPPLER  DATA 
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A single  radar  system  can  only  measure  the  radial  velocity  of  the 
scatterers  relative  to  the  radar.  This  component  of  the  scatterer 
motion  is  not  sufficient  to  characterize  the  three-dimensional  motion  of 
the  scatterers  within  the  sampling  volume  defined  by  a range  resolution 
element  times  the  antenna  beam  cross  section.  Models  must  be- employed 
to  extract  useful  data  from  the  Doppler  velocity  estimates.  If  a radar 
were  completely  surrounded  by  scatterers  all  moving  in  the  same  direction, 
the  particle  velocity  could  be  measured  by  making  observations  in  three 
different  directions  (including  vertical  for  three-dimensional  motion) . 
Unfortunately,  the  scatterers  are  not  all  moving  in  the  same  direction 
especially  in  the  vicinity  of  the  small  active  cells. 

4.1  Velocity  Information  Associated  with  a Small  Cell 

The  flow  field  about  and  within  a small  cell  is  too  complex  to  be 
measured  with  a single  Doppler  radar.  Shear  values  can  be  calculated 
for  the  area  within  the  cell  to  characterize  the  variation  of  the  flow 
field  within  that  cell.  For  a simple  axisymmetric  flow  pattern  model 
with  a vertical  symmetry  axis  radial  shear  values  can  be  identified  with 
convergence  and  tangential  shear  values  with  vorticity  if  measurements 
are  made  at  a low  elevation  angle.  The  success  of  the  plan  shear  indi- 
cator (Donaldson,  1970)  and  the  use  of  the  mesocyclone  signature  (Burgess, 
1976)  and  the  tornado  vortex  signature  (Brown  and  Lemon,  1976)  are  based 
upon  this  model  for  the  flow  field.  In  general,  the  flow  pattern  is  not 
axisymmetric  and  larger  scale  shear  deforms  the  simple  model  flow  causing 
a more  complex  pattern.  The  deviations  from  the  model  appear  to  be  small 
and  the  model  seems  to  be  useful  for  identifying  potential  regions  for 
the  development  of  tornadoes.  For  this  application  the  average  shear 
values  within  a cell  are  of  interest. 

Detailed  reflectivity  and  Doppler  velocity  measurements  in  situations 
characterized  by  supercells  and  tornadoes  reveal  mesoscale  circulation 
patterns  not  within  the  confines  of  a small  cell.  (See  the  Stillwater 
tornado  data  reported  by  Zrnic  et  al,  1976  and  discussed  in  Section  6.1; 
see  also  Agee  et  al , 1976.)  The  circulation  about  a weak  echo  region 
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appears  to  be  associated  with  a secondary  circulation  caused  by  the  cells  • 

neighboring  the  weak  echo  or  echo-free  vault.  The  mechanism  for  triggering  i 

and  maintaining  the  circulation  is  uncertain.  It  is  evident  that  a tangen- 
tial shear  signature  occurs  that  is  not  within  a cell.  This  region  can  | 

be  separately  identified  using  Doppler  velocity  data.  j 

f 

4.2  Mean  Velocity  Within  a Larger  Echo  Region 

The  lower  reflectivity  regions  surrounding  the  small,  active  cells 
generally  follow  the  environmental  wind.  Velocity  measurements  made  in 
the  lower  reflectivity  regions  may  be  used  to  estimate  the  environmental 
wind.  Observations  must  be  made  at  the  same  range  and  at  least  at  two 
different  azimuth  angles.  An  estimate  can  be  generated  by  assuming  that 
the  wind  at  a given  height  (range)  is  constant  over  the  azimuth  span  of 
each  echo  region.  The  velocity  is  calculated  using  a least  squares 
procedure  on  all  data  not  included  in  detected  cells  and  in  regions  with 
high  velocity  variance. 

4.3  Turbulence  Estimates 

Analysis  of  aircraft  observations  of  wind  velocity  fluctuations 
within  thundershowers  show  that  the  turbulence  is  anisotropic  at  scale 
sizes  in  excess  of  200  meters  and  suggest  that  velocity  variance  measure- 
ments at  scale  sizes  larger  than  200  m will  not  describe  the  turbulent 
dissipation  process  (Crane,  1976).  The  doppler  velocity  fluctuations  are 
caused  primarily  by  radial  velocity  fluctuations  at  scale  sizes  the  order 
of  the  antenna  beam  cross  section  at  the  measurement  range.  For  most 
radar  systems,  the  scale  sizes  associated  with  Doppler  measurements  are 
in  the  1 to  3 km  range,  significantly  outside  the  range  for  isotropic 
turbulence.  The  variance  estimates  therefore  correspond  to  larger  scale 
processes  such  as  organized  up  and  downdrafts  and  their  associated  con- 
vergence and  rotation  patterns,  for  example,  mesocyclones . The  radial 
wind  speed  profile  however  does  not  vary  linearly  across  the  radar  beam 
and  simple  models  to  estimate  variance  due  to  shear  will  lead  to  large 
measurement  errors.  Errors  in  a simple  linear  model  or  any  other  model 
for  the  variation  of  wind  speed  across  the  beam  will  cause  insufficient 
estimation  accuracy  to  remove  the  effect  of  shear  and  leave  an  accurately 
estimated  residual  component. 
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Pulse-volume  to  pulse-volume  changes  in  the  mean  Doppler  velocity 
will  be  used  to  estimate  shear  because  the  Doppler  variance  estimates 
should  be  identified  with  shear  and  because  the  variance  estimates  tend 
to  be  biased  and  severely  affected  by  noise.  The  radial  velocity  varia- 
tions within  the  beam  that  contribute  to  the  observed  variance  may  be 
associated  with  either  vertical  or  horizontal  variations  in  the  wind 
field.  Mean  Doppler  velocity  observations  may  be  used  to  estimate  the 
horizontal  variations  evident  after  averaging  by  the  antenna  beam.  Com- 
parisons should  be  made  between  tangential  shear  and  velocity  variance 
data  to  determine  the  degree  to  which  the  larger  scale  horizontal  shear 
contributes  to  the  variance.  If  significant  variance  can  occur  in  regions 
where  the  tangential  shear  is  low,  relative  variance  maxima  attributes 
could  also  be  used  to  characterize  the  radar  data. 
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5.  SOFTWARE  DESCRIPTION 


5 . I Program  Structure 

The  object  of  this  contract  was  to  develop  a set  of  algorithms  for 
the  processing  of  single  Doppler  radar  data  and  to  provide  a computer 
program  to  accomplish  the  data  processing.  The  cell  detection  procedure 
selected  for  this  task  is  based  upon  a procedure  previously  developed  by 
Crane  (1976).  The  algorithms  developed  under  this  contract  are  signifi- 
cantly different  from  the  earlier  ones  used  by  Crane  or  from  the  contour- 
ing algorithms  generally  used  on  large  scale  computers.  The  new  algorithms 
were  developed  specifically  for  this  contract  to  provide  rapid  computer 
processing  requiring  a minimum  of  computer  storage.  The  algorithms  were 
also  generated  to  simultaneously  process  both  reflectivity  and  Doppler 
data  in  a manner  constant  with  ultimate  employment  in  real-time  programs 
on  a mini -computer  coupled  to  a weather  radar. 

The  computer  program  processes  the  digital  radar  data  and  generates 
the  cell  and  larger  echo  area  attributes  identified  by  asterisks  in  Table 
1.  The  program  was  designed  to  read  digital  radar  data  tapes  prepared  by 
the  Weather  Radar  Branch  of  AFGL  at  their  Sudbury  field  station.  The 
raw  data  consisted  of  received  power,  mean  radial  velocity,  and  velocity 
variance  values  together  with  radar  operating  and  pointing  parameters 
plus  time  as  described  by  the  input  data  format  given  in  Appendix  A.  A 
series  of  subroutines  were  developed  to  read  and  reformat  the  radar  data, 
find  both  fixed  and  peak  referenced  contours,  and  calculate  the  attributes 
associated  with  the  contours.  A schematic  of  the  program  is  given  in 
Figure  4.  The  program  provides  contour  output  data  for  input  to  a 
second  program  which  generates  plots  for  fixed  level  contours,  and  out- 
puts attributes  calculated  for  the  small  convective  cells,  tangential 
shear  maxima,  and  fixed  echo  regions. 

This  program  is  configured  to  be  the  first  in  a series  of  programs 
that  (1)  detect  the  cells  and  generate  the  lists  of  attributes;  (2)  com- 
bine data  from  separate  scans  within  a volume  scan  to  provide  the  verti- 
cal development  attributes  for  the  detected  cells,  tangential  shear  maxi- 
ma, and  fixed  contours;  and  (3)  combine  data  from  separate  volume  scans 
to  generate  cell  tracks  and  to  list  the  time  histories  of  the  cells.  A 
schematic  overview  of  the  entire  processing  sequence  is  given  in  Figure  5. 
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The  computer  program  listing  is  reproduced  in  Appendix  B;  the  operating 
instructions  are  in  Appendix  A. 


5.2  Contour  Generation  Algorithm 

The  contouring  algorithm  used  to  find  both  the  fixed  level  contours 
and  the  peak  referenced  contours  was  designed  to  process  the  radar  data 
a single  radial  (all  data  for  a single  pointing  angle)  at  a time.  The 
processing  algorithm  was  tailored  after  the  technique  generally  used  to 
obtain  isoecho  contours  for  a weather  radar  display  and  is  significantly 
different  from  the  edge  following  algorithms  generally  used  in  computer 
processing.  The  edge  following  contouring  algorithm  requires  the  stor- 
age of  the  entire  data  field  in  the  main  computer  memory  at  one  time. 

For  the  radar  data  to  be  processed,  the  reflectivity  data  alone  would 
require  184,320  storage  locations  which  exceeds  the  available  core  stor- 
age if  not  packed  into  the  CDC  6600  computer  words.  If  the  data  were 
packed,  considerable  time  would  be  expended  unpacking  the  data  for  use 
with  the  contouring  algorithm.  This  new  approach  was  taken  to  minimize 
both  the  computer  storage  and  time  requirements.  The  processing  is  per- 
formed in  the  range,  azimuth  coordinates  of  the  radar.  The  program 
searches  the  data  in  range  along  a radial  defining  regions  or  events 
where  the  data  exceed  the  thresholds  for  contouring.  For  fixed  level 
contouring,  the  thresholds  are  preselected;  for  peak  reference  contour- 
ing, the  thresholds  are  computed  from  the  data.  This  contouring  algorithm 
differs  from  the  usual  application  of  isoecho  contouring  techniques  by 
combining  data  for  each  event  from  one  radial  to  the  next  to  generate  the 
attributes.  The  peak  detection  algorithm  is  unique  since  it  stores  suf- 
ficient data  from  radial -to-radial  to  obtain  the  required  attributes 
even  though  the  threshold  level  is  not  known  apriori. 

The  contouring  operation  starts  by  searching  the  data  along  a single 
radial.  The  start  and  stop  ranges  for  each  event  are  defined  by  threshold 
crossings  as  illustrated  in  Figure  6.  The  data  are  quantized  prior  to 
contouring  and  the  thresholds  are  applied  just  above  the  reported  value. 

For  example,  a 20  dBZ  threshold  would  include  only  data  that  exceeded  f] 

20  dBZ.  Since  a round-up  operation  is  included  in  the  generation  of  the  I 

quantized  data,  a 20  dBZ  threshold  would  include  all  values  above  20.5  I 

dBZ.  The  data  are  searched  by  threshold  at  each  range  element  reducing  I 
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the  number  of  tests  applied  at  each  range  element  to  a minimum.  The  ] 

event  identification  algorithm  is  depicted  by  the  flow  chart  on  Figure  7. 

In  the  remainder  of  the  processing,  only  data  within  an  event  are  tested 
or  combined  with  other  data  to  generate  attributes. 

The  data  for  events  from  one  radial  are  combined  with  data  from 
events  for  the  previous  radial  to  calculate  the  attributes.  This  process 
is  illustrated  schematically  on  Figure  8.  Events  on  both  radials  are 
searched  to  locate  adjacent  events.  If  more  than  one  B event  (previous 
radial)  overlaps  a single  C-event  (current  or  this  radial)  then  the  attri- 
butes for  both  B-events  are  combined  into  a single  set.  Each  identifiable 
echo  region  is  tagged  by  an  identification  number  which  is  used  to  index 
the  final  set  of  attributes. 

The  attributes  are  processed  separately  for  each  threshold.  Addi- 
tional processing  is  performed  for  the  lowest  level  threshold.  Each 
separate  identifiable  peak  along  each  radial  is  located  and  recorded  for 
subsequent  use  in  the  peak  reference  contouring  subroutine.  The  height 
of  each  range  element  within  each  lowest  threshold  event  is  calculated 
and  then  used  to  index  arrays  for  accumulating  reflectivity  and  velocity 
data  as  a function  of  height.  These  data  are  used  for  the  generation  of 
reflectivity  and  environmental  wind  profiles. 

The  peak  reference  contour  algorithms  are  identical  with  those 
described  above  with  the  exception  that  the  contouring  thresholds  are 
separately  calculated  for  each  radial.  The  peak  detection  algorithm 
uses  a threshold  a fixed  number  of  quantization  steps  below  the  peak 
value.  Since  the  peak  value  is  not  known  apriori,  attributes  must  be 
summed  for  each  possible  cell  (segment  of  radial)  within  the  fixed  num- 
ber of  steps  below  each  peak  value.  The  data  are  processed  one  event 
(lowest  fixed  level  threshold)  at  a time.  Threshold  levels  are  esta- 
blished for  each  peak  within  the  event.  The  data  at  each  threshold  level 
are  associated  between  the  B-  and  C-radial  segments.  Cells  are  declared 
when  cells  have  been  detected  which  do  not  enclose  other  cells  at  a 
threshold  level  the  required  fixed  number  of  steps  below  the  peak  value. 

When  B-  and  C-radial  data  are  associated,  the  highest  peak  from  either 
B-  or  C-  is  taken  as  the  new  peak  and  athe  attributes  are  restored  so 
only  data  for  the  required  fixed  number  of  steps  below  each  peak  are  saved. 

This  process  is  repeated  from  one  radial  to  the  next  until  a cell  is  not 
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updated  and  no  higher  level  data  are  present  on  the  next  radial  adjacent 
to  the  cell.  At  this  point,  a peak  referenced  cell  has  been  detected. 

Only  the  attributes  for  the  lowest  saved  threshold  relative  to  the  peak 
are  then  saved  for  subsequent  processing.  To  ensure  that  a second  cell 
immediately  adjacent  to  a previously  detected  cell  is  not  subsequently 
detected,  the  C-radial  segments  are  also  compared  with  B-radial  data  and 
attributes  for  a threshold  are  saved  only  when  the  C-radial  data  are  of 
higher  value  or  a B-radial  cell  is  being  activity  processed.  The  peak 
detection  process  is  illustrated  schematically  in  Figure  9. 

5.3  Attributes 

The  area,  average  reflectivity,  and  centroid  location  are  calculated 
for  each  fixed  contour  echo  region  and  for  the  contour  a prescribed  num- 
ber of  quantization  units  (CDB)  below  each  peak  value.  The  basic  data 
were  obtained  in  a polar  coordinate  system.  The  attributes  are  calculated 
as  follows  when  the  sums  are  taken  over  all  i (range),  and  j (azimuth) 
enclosed  within  the  contoured  region: 

A =.J:.(9.  - 0.  Jr.Ar 
i,J  J J-1  1 

Z”  = X •^•(9-  - 0-  Jr.Z.  .At 
A 1,3'-  3 3-1''  1 13 


X = .Z.  (9.  - 0.  ,)r.^sin9.Z. .Ar 

AZ  1,3  '■3  3-1''  1 3 13 


y = .E.  (0.  - 0.  , )r. ^cos9 .Z. .Ar 

y AZ  1,3  ''  3 3-1''  1 3 13 


where  A is  the  area,  Z is  the  average  of  the  logarithm  of  reflectivity, 

X,  y are  the  rectangular  coordinates  of  the  centroid,  9^  is  the  azimuth 
angle,  r^^  is  range,  Ar  is  the  range  interval,  and  Z^^  is  the  logarithm  of 
the  reflectivity  value  (in  dBZ) . For  the  detection  of  localized  tangen- 
tial shear  maxima,  the  logarithm  of  the  reflectivity  value  is  replaced  by 


VS.  . = (V.  . - V.  . J/(r.  (0.  - 0.  J) 
13  13  13-I'"  ^ i*-  3 3-1'''' 
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where  is  the  tangential  shear  and  is  the  mean  radial  velocity. 

Additional  shear  attributes  are  calculated  for  the  peak  reflectivity 
referenced  attributes.  These  were  the  average  radial  shear 

VR  = -J-.S.  (0.  - e.  ,)r.(V..  - V..  ) 

A 1.1  1 1-1''  i''  ij  ij-1'' 

the  average  tangential  shear 


and  the  average  radial  velocity 


V = 


i E. 

A 1,1 


(0.  - 9.  Jr.V.  .Ar 
1 1-1  1 11 


The  fixed  contour  profiles  are  calculated  by  summing  the  required 
attributes  for  height  regions  quantized  in  1 kilometer  steps.  The  height 
is  computed  using 
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r.sina  + 
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.2  2 
rj  cos  g 
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where  a is  the  elevation  angle,  R the  radius  of  the  earth  and  a the  effec- 
tive earth's  radius  multiplier.  A value  of  1.21  was  used  for  a.  In  addi- 
tion the  environmental  wind  velocity  profile  is  statistically  calculated 
using  mean  radial  velocity  data  confined  to  a narrow  reflectivity  interval 
(typically  20  to  35  dBZ)  and  for  sampling  elements  with  velocity  variance 
values  below  a preset  threshold,  TS.  The  mean  easterly  and  northerly 
velocities  u,  v are  calculated  as  follows: 

= [ (. Z .cos^0 .) (. E .V. . sin0 . ) - (. Z . sin9 . cos0 . ) (ZV. .cos9 . ) l/DEL 
1.1  3 1,3  13  3 1.3  1 1 13  1 

v(H)  = [(  Z.  sin^0.)  (ZV  cos0.)  - ( X sin0  . cos0  .)  (ZV.  . sin0  .)  1/DEL 
13  1 ij  3 13  3 3 13  13  3 

DEL  = (X.  sin^0  .)  ( J. cos^0 . 1 - (X. sin0  .cos0 
13  3 13  3 13  3 3 


where  is  the  mean  radial  velocity  and  the  summations  were  again  taken 
only  over  the  area  within  an  event  (identifiable  larger  echo  region). 


5.4  Sample  Results 


Processing  for  the  C-band  Doppler  radar  at  the  AFGL  Weather  Radar 
Branch  in  Sudbury  was  accomplished  using  512  range  intervals  of  300  m 


each.  The  raw  data  were  averaged  to  reduce  the  original  1024  range  ele- 
ments to  the  final  512  value.  The  processing  program  is  flexible  in 
adjusting  to  the  angular  increment  between  radial s.  For  the  data  from 
Sudbury,  the  interval  is  roughly  1“.  If  the  entire  360  by  512  data  array 
were  stored  for  both  reflectivity  and  mean  radial  velocity,  368,764 
words  of  core  storage  would  be  required,  roughly  four  times  the  106,000 
words  available  on  the  CDC-6600  computer  at  AFGL.  The  computer  program 
described  above  performs  the  required  contouring  and  attribute  generation 
operations  within  the  core  storage  available  on  the  computer  and  also 
provides  computer  generated  plots  of  the  fixed  level  contours  and  the 
centroid  locations  of  the  detected  cells. 

The  operation  of  the  fixed  contour  and  peak  detection  algorithms 
can  be  summarized  by  the  following  synthetic  example.  The  data  to  be 
contoured  are  given  in  Figure  10.  For  this  example,  the  threshold  for 
fixed  contouring  is  a value  of  0;  all  numbers  shown  are  within  the  fixed 
contour.  Table  2 depicts  the  values  for  the  start  and  stop  ranges  (I), 
the  event  number  and  the  echo  area  identifier  that  is  determined  after 
B-radial,  C-radial  association.  In  this  example,  all  the  data  are  for 
one  echo  area  or  region  although  as  many  as  3 events  are  detected  on  a 
single  radial.  The  number  of  peaks  and  their  locations  within  an  event 
are  also  listed.  Each  column  corresponds  to  an  array  in  the  program; 
their  function  is  explained  in  the  description  of  the  contouring  algorithm 
(Section  5.2). 

The  operation  of  the  peak  detection  algorithm  is  summarized  by  the 
entries  in  Table  3.  The  azimuth  and  event  values  are  the  same  as  for 
Table  2.  The  thresholds  generated  for  each  of  the  peaks  as  well  as  the 
segment  start  and  stop  locations  and  associations  as  possible  cells  and 
detected  cells  are  listed.  Note  that  the  start  locations  are  of  the 
range  element  preceding  the  threshold  crossing  as  are  the  stop  locations. 
The  cells  detected  by  the  algorithm  are  indicated  by  the  solid  lines  in 
Figure  10.  In  many  cases,  a zero  is  listed  in  the  possible  counter 
column.  In  these  cases,  no  cell  attribute  updating  takes  place.  A cell 
is  detected  when  a cell  is  not  updated  on  the  current  or  C-radial  and  no 
higher  adjacent  values  are  present  on  the  C-radial. 

Program  operation  to  date  has  been  to  debug  and  evaluate  the  operation 
of  the  new  algorithms  developed  for  the  fixed  contouring  and  peak  detection 
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operations.  The  program  has  been  exercised  using  both  a complex  ground 
^ clutter  pattern  that  severely  tests  the  multithreshold  program  logic 

^ ' required  for  peak  detection  and  actual  rain  data.  The  Doppler  measure- 

ments  - shear  values  and  average  radial  velocity  values  - are  all  zero 
for  the  ground  clutter  providing  a reasonable  technique  for  ground  clutter 
I suppression.  Sample  program  outputs  for  rain  are  depicted  in  Figures  11 

; through  16.  Figures  11  and  12  depict  the  B-Scan  printout  display  avail- 

, able  from  the  program  as  an  option.  The  data  are  averaged  in  range  and 

printed  for  each  azimuth.  The  solid  radial  lines  on  Figure  13  depict  the 
I start  and  stop  scan  boundaries.  Ground  clutter  is  evident  at  short 

1 ranges  at  many  azimuths  and  regions  of  rain  are  evident  to  the  west  and 

• north  of  the  radar.  Figures  13  and  14  depict  20  dBZ  reflectivity  con- 

tours. The  contours  in  Figure  14  were  obtained  from  the  second  EXPAND 
program  which  allows  variable  plotting  scales.  It  depicts  an  expanded 
view  of  the  region  to  the  south  and  west  of  the  radar.  The  fixed  contour 
identification  numbers  are  shown  together  with  dots  to  indicate  centroid 
location.  Finally,  the  attributes  are  listed  in  Figures  15  and  16. 


6.  SUMMARY  OF  RESULTS  AND  RECOMMENDATIONS 

6.1  Use  of  Attributes 

The  parameters  to  be  estimated  using  radar  volume  scans  are  listed 
in  Table  1.  The  attributes  can  be  objectively  obtained  from  reflectivity 
and  mean  velocity  estimates  from  a pulse-pair  processor.  The  cell 
detection  algorithm  was  tested  by  Crane  (1976)  using  a threshold  of  3 dB 
and  a precision  of  0.5  dB  (128  independent  samples).  A cursory  examination 
of  the  measurement  precision  problem  suggests  that  adequate  operation 
could  be  obtained  using  a precision  of  1 dB  (32  independent  samples) . 

This  problem  still  must  be  considered  in  some  detail  using  live  data. 

Sample  results  have  been  generated  using  pulse-pair  velocity  and 
reflectivity  data  provided  by  the  National  Severe  Storms  Laboratory 
(Doviak,  1976).  These  data  were  taken  with  a precision  of  2 dB  rms. 

Even  though  the  data  are  not  as  precise  as  desired,  the  expected  rela- 
tionships between  reflectivity  and  velocity  attributes  were  evident. 

Figures  17  through  19  depict  data  from  the  Stillwater  tornado  (Zrnic  et 
al,  1976).  The  reflectivity  data  on  Figure  17  have  been  simplified  by 
including  only  two  fixed  contours  and  identifying  the  locations  of  the 
small  cells.  The  reflectivity  data  as  provided  by  Doviak  contained  only 
5 dB  interval  contours  and  the  true  number  of  cells  could  exceed  the  num- 
ber displayed.  The  40  dBZ  contour  displays  a hooked  echo  although  the 
hook  pattern  could  not  be  discerned  at  any  other  contour  level.  The 
Stillwater  tornado  occurred  within  the  hook.  Figure  18  displays  the 
important  features  of  the  mean  velocity  pattern.  The  cell  locations  are 
again  displayed  on  this  figure.  The  highest  (positive)  and  lowest  (neg- 
ative) Doppler  velocity  contours  together  with  the  contour  midway  between 
the  two  are  displayed  for  each  shear  maxima.  The  tangential  shear  values 
are  listed  on  the  figure.  The  two  shear  maxima  between  the  0 and  10  km 
differential  X positions  straddled  the  20  dBZ  contour  and  are  caused  by 
noise  and  calculations  contaminated  by  using  data  from  regions  without 
scatterers.  These  shear  values  are  not  real  and  should  be  ignored.  The 
highest  shear  value,  0.03  s \ corresponds  to  the  Stillwater  tornado. 

This  value  occurs  in  a weak  echo  region  that  does  not  coincide  with  a 
cell  but  appears  between  two  cells  5 km  apart.  The  other  region  of  high 
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tangential  shear  is  associated  with  a gust  front  near  the  surface.  It  is 
interesting  that  a new  tornado  subsequently  formed  in  the  rather  weak 
echo  region  at  the  point  of  highest  tangential  shear  along  the  gust  front. 
Figure  19  displays  the  regions  with  rms  Doppler  velocity  fluctuations  in 
excess  of  10  m s ^ together  with  the  cells  and  high  tangential  shear 
regions.  In  this  figure  the  data  associated  with  the  edge  of  the  echo 
region  has  been  suppressed. 

The  data  displayed  in  Figures  17  to  19  show  the  important  details 
of  the  reflectivity  structure  and  the  associated  tangential  shear  field. 
Agee  et  al  (1976)  recently  reported  on  multiple  tornado  occurrences 
within  a single  mesoscale  cyclone  which  suggest  that  the  tornadoes  are 
more  closely  associated  with  the  reflectivity  maxima  that  move  within 
the  larger  scale  flow  field  surrounding  the  weak  echoregion  than  with  the 
weak  echo  region.  These  results  indicate  that  the  spatial  structure  of 
highly  localized  tangential  shear  maxima  and  possibly  associated  small 
convective  cells  is  an  important  characteristic  of  severe  storms  that 
spawn  tornadoes.  The  computer  program  developed  under  this  contract 
provides  a means  to  obtain  the  significant  information  from  the  radar 
and  display  it  in  a form  that  is  easy  to  interpret. 

6.2  Use  of  Program 

The  object  of  this  contract  was  to  develop  a computer  program  that 
would  significantly  reduce  the  amount  of  data  required  to  characterize 
a set  of  weather  radar  observations.  The  computer  program,  although 
designed  for  use  on  the  CDC-6600  computer  was  to  be  easily  transferred 
to  smaller,  dedicated  radar  site  computers.  The  program  that  was  devel- 
oped uses  algorithms  that  minimize  computer  storage  requirements.  The 
inital  program  uses  nearly  the  full  106,000  words  available  on  the  CDC- 
6600  computer  and  techniques  have  already  been  devised  to  significantly 
reduce  this  requirement  without  significantly  increasing  operation  time. 

As  currently  configured,  an  operational  version  of  the  program  can  be 
generated  that  uses  less  than  64,000  16-bit  words. 

To  date,  program  operation  has  only  been  to  debug  and  check  the 
program.  The  processing  includes  the  generation  of  a considerable  amount 
of  intermediate  output.  Processing  runs  have  taken  less  than  one  second 
of  CDC-6600  time  per  radial.  For  real  time  operation  with  an  onsite  com- 
puter, operate  times  of  less  than  one  third  to  one  sixth  this  value 
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are  required  and  should  be  achieved  with  available  computers.  Initial 
estimates  for  the  time  required  just  to  read  in,  calibrate,  and  store 
the  radar  data  for  a single  radial  range  between  .1  and  .2  seconds,  a 
significant  fraction  of  the  computer  time  required  for  all  the  process- 
ing. Preprocessing  of  the  data  to  provide  range  correction  and  scaling 
will  also  be  important  for  reducing  the  program  cycle  time  to  provide 
a real  time  capability. 

6.3  Recommendations 

The  program  developed  under  this  contract  is  a first  step  in  the 
generation  of  an  automatic  data  processing  system  for  single  station 
Doppler  radar  data.  The  program  has  been  subjected  only  to  preliminary 
analysis  to  ensure  that  the  computer  code  is  correct  and  the  program 
operates  as  designed.  Two  tasks  now  remain:  (1)  evaluate  the  operation 
of  the  program  with  a large  amount  of  weather  radar  data  and  (2)  gener- 
ate the  next  level  programs  to  track  the  cells. 
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Figure  3 Simplified  Display  of  the  Essential  Data  Contained  in  Figure 
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Figure  5 Overall  Processing  Scheme 
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Event  2:  Threshold  1 
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Figure  7 Event  Identification 
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Figure  10  Example  of  Detected  Peaks  (see  Tables  2 and  3) 
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Figure  lib  B SCAN  for  Reflectivity  (see  Figure  12  for  calibration) 
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Figure  12  B SCAN  Codes  for  Radial  Velocity  and  Reflectivity 


FTX'-D  ''ONlOf'  «TT'’I'’UT 


no-p  'i\-_ LPnfTTCiM  TITAl 


n 

THPrsMTin 
(n:»  ?) 

APrn 
(KM*  »?) 

iPFLPrTT'/T’’' 

(TP?) 

FAFT 
( K'') 

<Kt^) 

’7  53  IP 

pcrQi p 

1 

?n 

21 . P 

-32.  7 

-^7,  ■■ 

4.3P 

. pO 

? 

20 

<♦0.20 

24.  5 

- 42.  ° 

-54, 

3 0.25 

. 76 

3 

20 

.:n 

21.0 

-44.1 

-51 

. 1 4 

. 44 

4 

20 

21.0 

-4'.n 

-:i . 

. 1 4 

• ' 4 

Figure  15  Fixed  Contour  Attributes  for  Contours  Displayed  on  Figure  14 


PFAK 

OPTFCTFO  rp.LL  ATTT’IJTPP 

avr^An- 

A7  -:7A5P 

35  O'J 

2FPLPPTT\;iTV 

LCCA  r T' M 

Apfa  foot  j3;  TM 

PA  DIAL 

'AT-P'JTIAL 

'H5A’ 

PIPI AL 

'/PLO:  ITY 

TO 

1 

(fP’) 

?-  7 

(K*^) 

^ , T 3,  C; 

(k'M) 

-5e. 

(w/S/KM) 

R -.00 

(.3/ I/KM) 
0.00 

(M/S) 

-3.73 

2 

3 

27.3 

21-  . 0 

3.  ’ -47,  t 

1.3  -3«.4 

-P3. 
-4  8. 

7 .02 

1 0.00 

0.0  0 
0.00 

-S.  43 
-3.00 

Figure 

16  Small 

Cell  Attributes 

for 

Contours 

Displayed 

on  Figure  14 

49 


Y Differential  Distance  (km) 


X Differential  Distance  (km) 

(0,0)  is  at  (24,  82)  km  from  radar 


Figure  17  Reflectivity  Structure  for  the  Stillwater  Tornado  at  1.5  km 
Height 
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Figure  19  RMS  Doppler  Velocity  Fluctuation  Structure  for  the  Stillwater 
Tornado  at  1.5  km  Height 


APPENDIX  A 


PROGRAM  OPERATION 

A. I Description  of  Input  and  Output 

Program  input  and  output  are  depicted  in  Figure  A1 . The  tape  input 
format  is  given  in  Table  AI , The  control  cards  are  discussed  in  section 
A2.  The  program  produces  (a)  tapes  of  computed  attributes  for  input  to 
a second  program  for  computing  volume  scans;  (b)  a plot  tape  is  generated 
that  can  be  stored  for  input  to  another  program  "EXPAND"  which  is  a gen- 
eral purpose  plotting  package  for  plotting  the  fixed  contours,  centroids, 
cell  identification  and  peak  locations  expanded  over  selected  areas;  (c) 
B-scan  maps  are  also  produced  as  an  option;  (d)  full  scan  fixed  contour 
plots  of  the  lowest  threshold  level  can  be  obtained  on  35  mm  film  as  the 
program  is  executing;  (e)  hard  copy  plots  can  also  be  obtained;  and  (f) 
at  the  completion  of  a scan  the  program  will  print  out  fixed  contour 
attributes,  peak  detected  cell  attributes  and  tangential  shear  maxima 
attributes.  All  of  the  attributes  printed  have  identifiers  which  can  be 
associated  with  the  identifiers  displayed  on  the  expanded  plots. 

A. 2 Control  Card  Format 

Control  card  input  to  the  program  is  NAMELIST  input  which  allows 
certain  parameters  in  the  program  to  default  or  to  be  set  to  different 
values.  The  variable  names,  type  (LOGICAL  L,  INTEGER  i,  and  REAL  R) , 
dimension,  default  value  and  their  meanings  are  listed  in  Table  A2. 


TABLE  A2 


CARD  FORMAT  FOR  PROGRAM  EXTRAD 


Reads  in  program  parameters  via  NAMELIST  format. 


NAMELIST 

VARIABLES:  (Level 

760916) 

NAME 

TYPE 

DIMENSION 

DEFAULT 

MEANING 

PRINT  1 

L 

1 

FALSE 

When. True. Program  print  outs 
unpacked  raw  digital  data  from 
the  Doppler  data  tape. 

PRINT  2 

L 

1 

FALSE 

Currently  unused. 

PRINT  3 

L 

1 

FALSE 

When.True.B-Scan  maps  are 
produced. 

PRINT  4 

L 

1 

FALSE 

When. True. Full  scan  plots  are 
generated. 

ICODES 

I 

36 

A thru  Z then 

1 thru  9 fol- 
lowed by  a dot. 

Codes  for  representing  DBZ 
categories  for  B-Scan  map 
output . 

A1 

R 

1 ■ 

.13779 

In  the  linear  equation  y = mx+b 
For  computing  coded  DBZ  for 
B-scans,  A1  = M and  B1  = b. 

AZ 

R 

1 

.017 

Not  currently  used. 

BZ 

R 

1 

18.6 

Not  currently  used. 

CONTRZ 

L 

1 

FALSE 

IVhen. True. Fixed  contours  are 
generated  and  their  attributes. 
.False. will  ignore  fixed  con- 
touring . 

CONTRV 

L 

1 

FALSE 

When. True. Peak  detection  and 
their  attributes  will  be  gen- 
erated. .False. will  ignor  peak 
detection. 

NFILE 

I 

1 

1 

Not  currently  used. 

NUMF 

I 

1 

1 

Not  currently  used. 

AC 

R 

4 

-107. 7, +1.97, 
•.094, +.0018 

Calibration  coefficients  for 
computing  DBM  below  a threshold 
XCUT.  (See  XCUT.) 

r 


r 


NAME 

TYPE 

DIMENSION 

DEFAULT 

MEANING 

CALM 

R 

1 

.332 

In  the  calibration  equation 
y = mx+b,  CALM  = M and  CALB  » b. 

XCUT 

R 

1 

10.0 

Threshold  value  that  determines 
which  equation  to  use  for 
calibration,  (linear  or  non- 
linear) 

CK 

R 

1 

10.0 

In  the  equation  for  computing 
DBZ,  hence  K+P+ZOALOGIO(S (I- , 5) ) 
•CL+.5  (1)  K = CK. 

ZMAX 

R 

1 

0.0 

Not  currently  used. 

VMAX 

R 

1 

0.0 

Not  currently  used. 

NREC 

I 

1 

1 

Not  currently  used. 

NUMR 

I 

1 

999 

Number  of  radials  to  be  pro- 
cessed. Use  default  value  when 
doing  full  scan. 

IRUN 

I 

1 

0 

Run  number  chosen  by  user. 

INC 

I 

1 

0 

Not  currently  used. 

TL 

I 

4 

20,30,40,50 

DBZ  fixed  contouring  thresholds. 

LT 

I 

1 

4 

Number  of  fixed  contour  thres- 
holds to  produce.  (0<LT£4) 

TDW 

R 

1 

0.0 

Not  currently  used. 

DN 

R 

1 

0.0 

Not  currently  used. 

STARTR 

I 

1 

1 

Not  currently  used. 

DELTR 

R 

1 

1.0 

Not  currently  used. 

INPRF 

I 

1 

3333 

Value  of  PRF  (Pulse  Repetition 
Frequency) . To  be  used  when  PRF 
cannot  be  obtained  from  the 
data  tape. 

SCALE 

R 

1 

1.0 

Scale  factor  for  drawing  fixed 
contours . 

AE 

R 

1 

1.21 

Constant  for  computing  heights 
of  cells. 

AA 

R 

1 

300 

Constant  for  computing  heights 
of  cells. 

NAME 

TYPE 

DIMENSION 

DEFAULT 

MEANING 

BB 

R 

1 

l.S 

Constant  for  computing  heights 
of  cells. 

XI 

R 

i 

\ 

0.0 

Frame  size  coordinates  for 
fixed  contour  plotting.  Less 
than  or  equal  to  8 inches. 

X2 

R 

1 

8.0 

Same  as  above. 

Y1 

R 

1 

0.0 

Same  as  above. 

Y2 

R 

1 

8.0 

Same  as  above. 

TV 

I 

1 

35 

Velocity  attributes  are  not 
computed  for  DBZ  greater  than 
this  value. 

TSV 

R 

1 

10^ 

Not  currently  used. 

LDV 

I 

1 

3 

Cell  detection  treshold  for 
reflectance  peaks. 

LTV 

I 

1 

2 

Cell  detection  threshold  for 
velocity  peaks. 
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COMPUTER  PROGRAM  LISTING 
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i 
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m 

C 

c 

c 

PROGRAM  EXTRAD(INPtlT,'0lITPUTiTAPE5«INPUT»TAPE6»0llTPUTi  TAPE  IPO, 

• TAPE2,0PTp0,0EBUGp0UTPUT) 

program  EXTRAD  BRT  no,  162 

VERSION  2,0  level  761119 

MAIN  PROGRAM  SECTION, 

JHW  APGL  COC  6600 

logical  PR INT 1,PR1NT2, prints, PRlNTfl,CONTRZ,CONTRV 

integer  CmRD(S) 

COMMON  /PARM/  PRINTl,PRINT2,PRlNT3,PRlNTa,  ICnDE8(36)  ,A1,01,  A2,B2,' 
10nTRZ,C0nTRV,NPILE,NUMP,nREC,NUmR 

OATA  CwRD/4MPARA,aHEXEC,aHC0MM/ 

00000001 

00000002 

00000003 

00000004 

00000009 

00000006 

*00000007 

ooooooos 

00000009 

COOOOOOlO 

00000011 

00000012 

•00000011 

CALL  oat 

00000014 

1 

READ  C5,in  KEY 

00000019 

11 

PORMAT  (Att) 

00000016 

IF  (E0F(9))  91,21,91 

00000017 

21 

CALL  PAGE 

oooooois 

WRITE  16,31)  key 

00D00019 

31 

FORMAT  (IH  ,Afl) 

00000020 

00  41  KPl.B 

00000021 

IF  (KEY,EQ,CWRO(K))  GO  TO  (61,71,81),  K 

00000022 

ai 

CONTINUE 

00000023 

WRITE  (6,51) 

00000024 

91 

FORMAT  (16M  illegal  KEYWORD) 

00000029 

GO  TO  91 

00000026 

c 

00000027 

C 

* PARAMETERS  * PACKAGE. 

00000028 

C 

00000029 

61 

CALL  INPARM 

00000030 

GO  TO  1 

00000031 

C 

00000032 

C 

* execution  * package. 

00000033 

C 

00000034 

71 

CALL  EXTRAT 

00000039 

GO  TO  1 

00000036 

C 

00000037 

c 

* COMMENTS  CARD  * PACKAGE, 

00000036 

c 

00000039 

«1 

CALL  INE  (5) 

00000040 

GO  TO  1 

00000041 

c 

00000042 

c 

end  of  JOS, 

00000043 

c 

00000044 

•Jl 

write  (6,101) 

00000049 

101 

FORMAT  (//2X,7H  EnOJOB) 

00000046 

IF  (,N0T,PRINT4)  CO  TO  111 

00000047 

CALL  EnDPLT 

00000048 

111 

STOP 

00000049 

INO 

00000090 

o o o o 


BLHCK  OATA  00000051 

***«****«**********«******************ft*************************;*00000052 

FOR  PROGRAM  EXTRAO  ERT  lfc2  00000053 

VERSION  2,0  level  761119  0000005« 

JHW  C0C660O  AF6L  00000055 

******************************************************* ***********00000056 
logical  PRINTl , print?, PR1NT5;pRINT6,C0nTRZ»C0NTRV  00000057 

integer  TL»8TARTH,TV,T8V  00000058 

c .—.......•........•.•............................•00000059 

COMMON  /PaRM/  PRINT1,PRINT2,PRTnTS,PRInT4,  IC0DEa(36}  »A1,B1,  A2,B2.'C00000060 
10NTRZ»C0NtRv*NFILE.NUMF,NREC,NUMR  00000061 

common  /INSUB/  TL(«) ,LT,T0w,0N,8TARTR,0ELTR,RN(a) ,8C0N,CElnTH{5)  0000006? 
COMMON  /AZM/  AZMUTH(a60) ,NA,ELEvAT,PRF;kEEP  00000063 

COMMON  /AlO?a/  MVP(3,10?«)  00000064 

common  /VALMAX/  ZMAX,VMAX,ACf4) , calm, CALB,XCUT,CK, INC  00000065 

COMMON  /ADATA/  IOAV, IMOUR, IMIN, I8EC,NTP,N8F,n00,NRC  00000066 

COMMON  /HEAD/  TlTLE(6),IC0nE,VER8, level, DAT, IRON, NPAGE,NL0G  00000067 

COMMON  /LINUM/  line  00000068 

c ******************************************************************00000069 

COMMON  /MOPED/  INPRF, scale. LOV, LTV  00000070 

common  /store/  AE, AA,BB,SLtCL»Tv,T8V  00000071 

COMMON  /EXPAN/  XI,X?,Y1,Y2,XMIN,XMAX,YmIN,YMAX  00000072 

c ..................................................................00000073 

DATA  PR1NT1/,FALSE,/,PRIN72/,FAL8E,/,PRINT3/.FALSE,/»PRINt4/,FAL8E00000074 
1 ./,  Al/,13779/,Bl/l,5/,  A?/.017/,B?/18,6/,C0NTRZ/,FALSE,/,C0NTRV/,FA00000075 
2L8E,/»NFIlB/1/,NUmF/1/,N9EC/1/,NUMR/999/  00000076 

DATA  TL/20,30,40,50/  00000077 

DATA  X1/0,0/,X?/8,0/;y1/0,0/, Y?/8,0/,AE/l,2l/, AA/300,/,8B/1,5/  00000078 

DATA  3L/0,0/,CL/0,0/,TV/33/,T8V/1000000/  00000079 

DATA  LT/4/,8TARTR/1/,DELTR/1,0/  00000080 

DATA  TDW/0,0/,DN/0,0/,PN/256.0,512,0,768,0, 1 024,0/, scon/299, 79?5/  00000081 
DATA  CELMTM/0, 5, 1,042,2,0/  00000082 

DATA  IC00E8/1HA, IHR, IHC, IHO, 1 HE , 1 HF , 1 MG, 1 HH, 1 H 1 , 1 H J , IHK, IHL, IHM, 1H00000083 
IN,  IHO,  IHP,  IHO,  IHR,  IMS,  IHT,  IHIJ,  IHV,  IHW,  IHX,  IHY»  IHZ,  IHl , 1 M2 , 1 H3,  1 H4 , 00000084 
21H5, 1M6, IH7, 1H8, 1H9, IH,/  00000085 

DATA  ZM  A X /O,  0/,  VM  AX /0,0/,  AC /•  1 07,76555,’ I ,97  678  38,  ••094297528,  ,000100000086 
l8226318/,CALM/0,332/,CALB/*98,3/,XCUT/10,0/  00000087 

DATA  TITLE/7HPR0GRAM,7H  EXTRA0,1M  ,1H  ,IH  ,1H  / 00000088 

DATA  I9UN/0/,NPaBE/0/,IC00E/16?/, vERS/1 , 0/ , LEVEL /760916/  00000089 

DATA  InPRF/3333/  00000090 

DATA  CK/10,0/  00000091 

DATA  SCAUE/1  ,0/,LDv/3/, LTV/2/  00000092 

C ..................................................................00000093 

end  00000094 
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r>  o o ni  ri  r>  r»  o r>  r»  r»  c* 


SUBROUTINE  INPARH 


VERSION  2,0 
JHN  AP6L 
CONTROL  CARO 


level  T6111R 
C0C6600 

INPUT  parameters. 


logical  PRINTI , print?, prints, PR I NTR.CONTRZiCnNTRV 
integer  TL,8TARTR,TV,T8V 


OOOOOORS: 

tooooooR*; 

0000009T: 
00000098 
00000099; 
[00 
iOi; 
[02 
OS 
[04 


lOOOOO 

00000 

00000 

.00000 

00000 


,NLr)G 

,CElWTH(3) 


OIMEnSiON  PROGIDtS) 

COMMON  /PARM/  PRINTl, PRINT?, prints, PRINT4,IC0DE8(3<iJ, A 

iontrz,contrv,nfile,'nump,nrec,numr 

COMMON  /VALMAX/  ZMAX , VM AV , AC (4) , C ALM , C ALB, XCdT, CK , INC 
COMMON  /HEAD/  TITLE(8) , IC0nE,VER8, level, OAT, IRUN,nPAGE 
COMMON  /INSUB/  TL(a) ,LT,T0w,0N,8TARTR,nELTR,RN(«) ,8C0N 
COMMON  /MOREO/  INPRF,8CALE,LDV,LTV 
COMMON  /EXPAN/  XI , X?, T I , T?, XMIN, XMAX, YMIN, YMAX 
COMMON  /STORE/  AE,AA,B8,SL,CL,TV,T8V 

OATA  PROGIO/THWILLANO, IH  ,1H  / 
namelist  /input/  PRINTl,PRINT2,PRINT3,PRlNTa, ICODES, A1 , Bl , A2, B?, COOOOOO 
INTRZ,C0NTRV,NFILE,NUMF, AC,CALM,CAL8,XCUT,CK,2MAX, VMAX,NREC,NUMR,IROOOOO 
2UN,INC,TL,LT,T0w,DN,STARTR,0ELTR,INPRF;sCALE,AE, AA,B8,XI,X2, Yl, Y?,00000 


l,Bl,A?,B2,COOOOO 
00000 
00000 
00000 
00000 
00000 
00000 
00000 
00000 
00000 


3TV,T8V,LDV,LTV 

00000 

READ  (5, INPUT) 

00000 

IF  CE0FC5))  111,1,111 

00000 

WRITE  (6, INPUT) 

00000 

IF  (.NOT.CONTRZ)  GO  TO  21 

00000 

IF  (,N0T,PRINT4)  CO  TO  11 

00000 

CALL  CRTPLT  (PROGin, 1,0, 17,0) 

00000 

CALL  FLOT  (0,0, 0,0, 3) 

00000 

CALL  PLOT  (8,,0,,2) 

00000 

CALL  FLOT  (8., a, ,2) 

00000 

CALL  PLOT  (0,,8,,2) 

00000 

CALL  FLOT  (0,,0,,2) 

00000 

XiSIN(0.0)«4.0 

00000 

viC08Co,n)*e,o 

00000 

CALL  FlOT(X,Y,3) 

00000 

Y«Y-,25 

00000 

CALL  FL0T(k,Y,2) 

00000 

8CALE«8,0/(y2»YI) 

00000 

IF  ((X2.Xl),CT,(Y2-Yl))  SC ALE-S . 0/ ( X2-y 1 ) 

00000 

*mIn«8CALE«X1 

00000 

XMAX»SCALE*X2 

00000 

YMIN«8CALE*Y1 

00000 

YMAX«SCALE*Y2 

00000 

IF  (, NOT, PRINT?)  GO  TO  61 

00000 

00000 

print  icooes  values. 

ooooo 

00000 

CALL  PAGE 

ooooo 

write  (6,31) 

ooooo 

11 


IHVALUE/40X,7HAN0  PWR) 


FOR  VAR, TX, 500000 
00000 


09 

06 

or 

06 

09 

10 
11 
12 

13 

14 
19 
16 
IT 
18 

19 

20 
21 
22 
29 
24 
29 
26 
2T 
28 

29 

30 

51 

52 
S9 
34 
39 
36 
ST 
SO 

39 

40 

41 

42 

43 

44 
49 
46 
4T 

48 

49 


63 


on  41  I»l,36 

00000180 

XA«(FL04T(I).B1)/AI 

00000151 

IF  CXA.lt, 0.)  XA»0. 

00000182 

X8«(FLa4T(I)*B2)/48 

OOOOOISI 

41 

write  (6.51)  ICQDE8(n/xB*lC00E8(n»XA 

00000184 

51 

FORMAT  C15X,Al,9J{,FR,3;ilX,Al,9x.F9,3) 

oooooiss 

&1 

CONTINUE 

0000018B 

IF  (.NOT, PRINTS)  GO  TO  lOl 

00000187 

CALL  PAGE 

00000188 

write  C6,Tl) 

00000189 

Tl 

FORMAT  (IH0,8X, 13HC0DE  FOR  OBZ, TX» 5HVALUE) 

00000160 

00  81  I«1,S6 

00000161 

XAnCFLOATCn-BD/Al 

00000162 

IF  (XA.LT.O.)  XABO, 

00000163 

SI 

write  (6»R1)  IC00E8(I)«XA 

00000164 

41 

FORMAT  (15X,Al,9X,F9,3) 

00000168 

101 

CONTINUE 

00000166 

RETURN 

00000167 

' 111 

write  (8,121) 

00000168 

121 

FORMAT  (30H  END  OF  FILE  IN  NAMELIST  INPUT) 

00000169 

STOP 

00000170 

End 

00000171 

o o o r>  o o 


SUBROUTINE  EXTRAT 


PJP  MODIFIED  4/27/77 
VERSION  1,0  LEVELP7ti0916 
JHW  AF6L  CDC  6S00 
UNPACKING  ROUTINE. 


MOO,  1,0 


LOGICAL  PRINT!, PRINT2 
integer  VJ(S14)  ,111(51 
INTEGER  W(514) , V(5ia) 
I T(50) ,TC(10,30) ,T0(1 

210.30) ,1P02(1O,IO,3O) 
3,lPCS(lO,10,30!,ICl(a 
4C(a,10,4),IPTC(30),TC 

50.10.30) »IPC3T(10,10, 
fe(10,SO) , IPRnG(10,30) , 

integer  TL,STARTR 
REAL  UPC8, 100) ,01(501 
i,4),ZH(12,15,100),n8l 
2VATR(17, 100) ,BI (3,50, 


00000172 

►00000173 

00000174 

00000175 

OOOOOlTfc 

00000177 

(00000178 


DATA  MVPriSK/77770000000000000000B, 
C000077770000000000008, 

C000000007 777000000008, 

C 00000000000077 77 00008, 
C00000000000000007777fi/ 

DATA  OAYM8K/00170000000000000000B, 

C 036000000000000000008, 

C7 40000000000000000008/ 

DATA  HRM8K/00000 1700000000000008, 

C 000006000000000000008/ 

DATA  M I NMSK/ 000000000 170000000008, 

C00000000160000000000B/ 

DATA  SECMSK/ 00000000000000 1700008, 

C 0000000000000 16000008/ 

DATA  lREGM8K/OnQO7777OOOOOOOOOOOO0/ 
DATA  N8H8K/00000000000000000600B/ 
DATA  VPH8K/00000000000000003770B/ 
DATA  ME ANM2/ 0000000000000000 17778/ 
DATA  I8ONMSK/OOOOOOOOOOOOOOOO2OOO0/ 
DATA  EUMSK/000000000000000077778/ 


,PR1NT3,PRINT4,C0NTRZ,C0NTRV  00000179 

4),vai(514)  00000180 

, vsr5l4) ,Sv(S14) ,VB(514),HB(514) ,HVB(514) , 00000181 

0,3O),KDD(4),IPT8(30),IPLn(10,l0),IPBl (10,  00000182 

,IP83( 10, 10,30) , IPCI (10,10,30) , IPC2(lO, 10,30)00000183 
,30,4),IPBNT(10,30) ,IPCNT(10,30) ,18(4,30,4) ,100000184 
1(10,10,30), IPTCl (30) ,IPC1T(10, 10,30), IPC2T (100000 1 85 
30) , ICVNT(4) ,I0VNT(4) , 10C(30) , I D VC (30) , IPVRNGOOOOO 1 86 
IPCNTT(10,30) ,IACT(100),IACV(100)  00000187 

00000188 

, HZ (2 , 15, 30) , VI (3, 15,30) » Cl (3, 30,4), atR (5, 10000000189 
(100) ,TATR(33,100) ,CI1 (3,30,4) ,UV(5, 100) , 00000190 

4)  00000191 

...........00000192 

00000193 
00000194 
00000195 
00000196 
00000197 
.'C00000198 
00000199 
00000200 
00000201 
00000202 
00000203 
00000204 
00000205 
.•00000206 
00000207 
00000208 
00000209 
00000210 
00000211 
00000212 
00000213 
00000214 
00000215 
00000216 
00000217 
00000218 
00000219 
00000220 
00000221 
00000222 
00000223 
00000224 
00000225 
00000226 


DIMENSION  IN(158) 

DIMENSION  OAYMSK  (3) ,IDAYSFT(3) ,HRM8K(2) , IHRSFT (2) , MINMSK (2) , 
CIMINSFT(2) ,SECMSK(2) ,18EC8FT(2) ,MVPM8K(5),IMVPSFT(5) 
COMMON/WORK/W, V, VS,8V, V8, VJ,UI, VSI,HB,HVB,NCL 
COMMON  /AZM/  AZMUTM(460) ,NA,ELEvAT,PRF,KEEP 

COMMON  /PaRM/  PRInT1,PRINT2,PRInT3,PRINT4, IC0DES(36) , A1,B1, A2,P2 
10NTRZ,C0NTRV,nFILE,NUMF,NREC,NUmR 
COMMON  /A1024/  MVP(3,1024) 

COMMON  /VALMAX/  ZMAX,VMAX,AC(4), calm, CALB,XCUT,CK, INC 
COMMON  /ADATA/  ID AV  , IHOljR , IM I N , 1 8EC , NTP , NSF  , NOD , NRC 
COMMON  /MOREO/  INPRF, scale, LDV, LTV 

COMMON  /INSUB/  TL(4) , LT , TDw , ON , ST ARTR, DELTR , RN (4) , SCON, CEL^TH (3) 
COMMON  /AZ2/8Ina,C0SA,DELTA2, I8CANF,NEL 
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iunon*«  *-  ooooooooo 


1 


I 

23 


25 


27 


data 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

lH/15/ 

DATA 


HE ANM 1/000000007777000000005/ 

V A RHSK 1/0000000000007 77 7 000 OB/ 

PwHHSKl/OOOOOOOOOOOOOOOOr  7 778/ 
AZMSK/0000777700000000COOOB/ 

NHSK/ 740000000000000000008/ 

»<RCH8K/0000000000000000003o8/ 

N8FMSK/00000000000000000600B/ 

NDOMSK/OOOOOOOOOOOOOOOOOOflOB/ 

NTPMSK/00000000000000003000B/ 

KM8K/00360000n0O00O<J000008/ 

iPARBSK/OOOOOOOOOOOOOOOnaOOOB/ 

IMVP8PT/12,2«,*24,*12,0/ 

I0AT8PT/12,8,a/ 

IMR8PT/21 , 17/ 
lHIMSFT/*27,29/ 

I8EC8ET/»12,*18/ 

JHAX/10/,KMAX/10/. IEMAX/30/, IAT/5/,NI0/l00/»NPC/4/,H7P/i 

,SUP/8/,N«Z/2/.NVI/3/#NPA/a/,NUHAX/33/,MTT/50/,NCL/5l4/ 

DAZT/2,0/,NP8/3/,NUV/5/,NVmAX/17/ 


lEMAX  ■ maximum  no,  EVCNTS/RAOIAL  t NId  ■ NO,  OF  IO'S/SCAn, 
wPa  ■ NO,  PARAMETER8  I NFC«NO,  FIXED  CnNT0UR8  I NA  ■ AZIMUTH 
LOB  • NO,  peak  contours, 

N0D«P»E0,  of  DUMP  PULSES  ALT»0»ALLm1 
NTP«CELL  NIDTH  0,1,2  MEANING  ,5,1,042,2, 
nSF»»iJBFRAmE  0, 1,2,3 

NRCBNO,  range  cells  0,1, 2, 3 MEANING  258,512,766,1 


00000227 
00000228 
00000229 
00000230 
00000231 
00000232 
00000233 
00000234 
00000235 
00000236 
00000237 
00000238 
00000239 
00000240 
00000241 
00000242 
2/,nZ0000024S 
00000244 
00000245 
....*00000246 
00000247 
00000248 
00000249 
00000250 
00000251 
00000252 
00000253 
00000254 


NO, 


024 


IEOFbO 

00000255 

ISCANFmO 

00000256 

NFC«LT 

00000257 

NAal 

00000258 

NEL»1 

00000259 

BUFFER  IN  (i,n  (lN(n,lN(15B)) 

00000260 

IF  (UNIT(I))  l,l8i;201 

00000261 

DO  11  I«l,5l4 

00000262 

W(I)»0 

00000263 

V (I) ■•999 

00000264 

V8(n»-999 

00000265 

SVf  n»*999 

00000266 

00000267 

UNPACK  DAT  HOUR  minute 

second  and  STATUS  FLAGS. 

00000268 

00000269 

IOAT«0 

00000270 

DO  23  !■!,  S 

00000271 

IDAT«I0AT^1  0«*(I»n  •SHIFT  CiNfl) 

,ANO,  OATMSK  f I) , lOATSFTd) ) 

00000272 

continue 

00000273 

IHOIJP^O 

00000274 

DO  25  I«l,2 

00000275 

lHnu9»IHOUR4lO**(I«nASHlFT(TN(i)  ,AND,  HRHSK  ( I ) , IHRSFT  ( 1)5 

00000276 

continue 

00000277 

IMINbO 

00000276 

do  27  Ib1,2 

00000279 

Imin»ImIN+10#*(I  IFTlINd) 

,ANO,  MlNM8K(n,lMlN8FT(I)) 

00000280 

continue 

00000281 
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o- r»  r»  o o o r»  r»ooo  uj « m o o o 


29 

I 


I 


I 

1 


[ 

f 

1 

I 


I 


isEC«o  oooooatt 
on  29  IM1,2  000002S1 
l8EC»l8EC*lO**(I.l)*8HIFTClN(n  .AND,  SECH8K 1 1 ) , I8EC8FT  ( I ) ) 0000028#; 
CONTINUE  0000028S 


NTF^SHIFTtlNd^  .AND,  NTPN8K,«9)  OOOOOE#*: 

NSF«8HIFT(IN(n  .AND.  NSMSK,«T)  0000028? 

nDD«8HTFT(1N(1)  .and,  N00N8K,-9)  00000288 

NHC«SHTFT(lN(n  .AND,  N>?CN8K,»1J  00000288' 

00000290 

UNPACK  PPF  azimuth,  and  ELEVATION,  00000291 


00000291 

00000291 

00000298 

0000029S 

00000298 

0000029? 

00000298 

00000299 

00000300 

OOOOOSOt 

OOOOOSOt 

00000303 

00000308 

00000309 

00000308 

0000030? 

00000308 

00000309 

00000310 

00000311 

OOOOOIlf 

0000031J 

00000318 


00000319 

UNPACK  THE  DATA, 

00000S18 

N8F  IS  SUBFRAME. 

0000031? 

00000318 

KbMSF*258*1 

00000319 

KEEPbK 

00000320 

00000321 

UNPACK  FIRST  DATA 

HORO, 

00000322 

00000323 

HVP(1,K)P8HIFT(IN(3)  .and. 

MEANMl,-24) 

00000328 

mVP(2,k)bSHIFT(IN(3)  .and. 

VARMSKl ,fl2) 

00000329 

MVP(3»K)i9HIFT(IN(1)  .AND, 

PWRMSKI.O) 

00000326 

N«3 

0000032? 

00000328 

UNPACK  REMAINING 

DATA, 

00000329 

00000330 

DO  63  IB1,1 

00000331 

NaN4l 

00000332 

DO  63  JBl,3 

00000333 

IF(I  .CQ.  1 .and.  (J  ,EQ. 

1 aOR.  J tEQ, 

4))KbK^1 

00000338 

IFCI  ,EQ.  1 .AND.  (J  .EG, 

1 .OR.  J .EO, 

4))MbI 

00000339 

IF(I  .EQ.  2 .AND,  (J  .EG. 

2 .OR.  J .EG. 

3))KbK*1 

00000338 

N»8WIFTCIN(2)  .AND,’  NM8K,tti 
K»9HIFT(IN(21  .and,  KM8K,in 
IF  (K.EO.O.AND.n.EQ.O)  go  TO  91 
JAbO 
J0bO 

DO  31  iBl.a 
lRE6«2**(l»n 

IF  CdREG.AND.N)  ,NF,0)  JA»JA+l 
IF  ((IREG.AnO.K) ,NE,0)  JPBja*! 

CONTINUE 

IF  (JA.EQ.l.ANO.JB.EQ.l)  PRFi394. 

IF  CJA.EO.l .AN0.J8,EQ,3)  PRFa794. 

IF  (JA,GE,3,AND.J8,LE,n  PRFbUH, 

IF  (JA,BE.3,AN0,J8,GE,3)  PRFb333S, 

GO  TO  31 
PRFbINPRF 

continue 

IRE6i8HlFT(IN(2)  .AND.  AZNSK,2a) 
AZMUTHtNA)BlREG*360,0/fl098, 
IREG*8HIFT(INC2)  .AND,  ELNSK.O) 

ELE V at* IREOb 380, 0/4096, 

IF  (BLEVAT.GT, 180,1  eL€VAT«ELEVAT-380, 
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IF(I  •CO,  2 .AND.  (J  .CO.  2 .OR,  J .CO.'  5))Mb1 

000003ST' 

IF(X  .CO.  3 .AND.  (J  ,10.  3))KiKB1 

000003SB 

IF(I  .EQ.  3 .AND.  (J  .CO.  S))M*i 

00000339 

MVP(M,K)b8H1FT(IN(N)  .and,  MVPM8K(J),IMVP8FT(J)) 

00000340 

M«M«1 

000003A1 

63 

CONTINUE 

00000341 

65 

continue 

00000343 

IF(N  ,LT.  155)  GO  TP  61 

00000344 

C 

00000345 

c 

clean  off  extra  bits. 

00000346 

c 

00000347 

00  71  iRKCEPfK 

00000346; 

IPARl»8HlFT(MVP(l,n  ,ANO,  IPARm8K,»11) 

00000349 

1REGbShIFT(HVP(1, I)  ,ANO,  I8CNMSK,-10) 

00000350 

MVP(1,X)RNVP(1,X)  .and,'  MEANM2 

00000351 

IF  (IREG,GT,0)  HVP(l,nB.NVP(l,I) 

00000352 

IPAR2»8HIFT(MVP(2, T)  ,ANO,  IPARm8K,*11) 

00000353 

mvP(1,I)i8H1FT(MVP(2,I)  ,AN0,  VPH8F,*3) 

00000354 

IPAR3R8HIFT(MVP(3,I)  ,AN0,  IPARhSK,*!!) 

00000355 

MVP(3,I)b8MIFT(HVP(3,I)  .AND.  VPHSK,ti3) 

00000356 

T1 

CONTINUE 

00000357 

IF  (•NOT.PRINTl)  60  TQ  81 

00000356 

CALL  PRNl 

00000359 

C 

00000360 

C 

GET  NEXT  TAPE  RECORD, 

00000361 

C 

00000362 

6i 

BUFFER  IN  (i,n  (IN(1),IN(158)) 

00000363 

IF  (UNITCD)  91,111,201 

00000364 

91 

nSb8HIFT(IN(1)  .and,  NSMSK,»7) 

00000365 

IREG«8HIFT(IN{2)  ,AND,  AZR8K,2«) 

00000366 

AzbIREG*360,0/4096. 

00000367 

DAZ«AZ»AZmutH(NA) 

00000366 

IF  (ABSCOAZ) ,6T,0AZT)  GO  TO  101 

00000369 

DAZSb8I6N(1,,DAZ) 

00000370 

IFfNA.EQ.nOAZFBOAZS 

00000371 

IF  (DAZF,CQ.OAZS)  00  TO  HI 

00000372 

101 

A2bA2b360,*DAZF 

00000373 

OAZaAZ»AZMUTH(NA) 

00000374  , 

IF  (ABSCDAZ) ,GT,DAZT1  SO  TO  121 

00000375  ; 

DAZS«8I6N(l.,0AZ) 

00000376  { 

IF  CDAZF,nE,0AZ8)  00  TO  121 

00000377  1 

IF  (A88(AZ«S60,*0AZF*AZMIJTH(1)),0T,DAZT)  00  TO  HI 

00000376  i 

C 

finished  scan 

00000379 

ISCAnFbI 

00000360 

AZRAZ«360.*0AZF 

00000361 

00  TO  131 

00000362 

111 

lEOF-l 

00000363 

121 

ISCANFb.I 

00000364 

131 

CALL  C0NTR2  ( JM AX , KM AX , lEH A X , I At, NIO , NFC , NZP, NZH, NUP, NUHAx , NHZ , NVI00000385 

1,NPA,IPvRnO,IPRnO,IOC,IOvC,IP81,IPB2,IP03,IPC1,IPC2,IPC3,IPTB 

,TB, 100000366 

2PBNT,IPCNT,T,IPTC,HP,TC,1B,IC,Mz,VX,ICvnT,IBvNT,DI,TCI,IPtCI, 

1PC17000003B7 

3,IPC2T,IPC3T,CI,ATR,ZH,D8I,IPL0,TATR,K0D,NTT,CIi,IC1,NPB, 

00000366 

«IPCNTT,IACT,IACV,N!IV,NVMAX,UV,  VATR,NNE,'BI) 

00000369 

IF  (IEOF,Ea,n  00  TO  181 

00000390 

IF  (NA,OT,NUMR)  RETURN 

68 

00000391 

1 


68 


lai 


NRCJCO  TO  21 


00000S92 
00000393 
00000394 

NZ2f NZH|NUP»NUMAy,NHZ|NVlOOOOOS9S 
, IPCI f lPC2i !PC3« IPTI.TB. 100000396 
, I0VNT,DI,TC1,IPtCI,1PC1TOOOOO39T 
NTT,CIl,ICl,NP8,  00000398 
) 00000399 


161 

NAbNA^I 

00000400 

IP  (NA,6T,»JUMR)  60  TO  121 

00000401 

GO  TO  1 

00000402 

181 

WRITf  (6,191) 

00000403 

191 

FORMAT  {19H  EOF  READ  ON  UNIT  1) 

00000404 

IF  (.N0T,PR1NT2)  60  TO  221 

00000405 

C 

CALL  PRN2{2) 

00000406 

60  TO  221 

0000040T 

201 

WRITE  (6,211) 

00000408 

211 

format  (21H  PARITY  ERR  ON  UNIT  1) 

00000409 

221 

RETURN 

00000410 

END 

00000411 

o o o MMonn  r»oo  rioo«-»  r»  ci  r>  o o r»  o o r»oor»o 


1 


1 

I 


subroutine  compz 


VERSION  1,0  LEVEL  TblllB 
JHW  AFGL  CDC6600 
CONFUTES  D8Z. 


integer  N, V, VS»SVf VB, VJfUl, VSIfHBfHVB 

integer  tl.startr 


C0MM0N/W0RK/W(5l«)  ,V(5ia),V8f5U)i8Vr5l4),VB(5l4),VJ(5ia)  ,UI(5ia) 
IV8I(5U)  »HB(514),HVB(514),NCL 
COMMON  /A1024/  MVP(3,1024) 

COMMON  /AZM/  AZMUTM(460)  ,NA,BLEv4T,PRF,’KEEP 

COMMON  /VALMAX/  ZMAX,VMAX»ACf4),CALH,CALB,XCUT,Ct<iINC 

common  /AOATA/  IDAY,IH0UR»IMIN,IBEC,NTP,N8F,n0D,NRC 

COMMON  /INSUB/  TL(4),LT,T0w,DN,STARTR,pELTR,RNf4),SC0N,CELWTM(3) 


FETCH  number  of  RANGE  CELLS  (N), 


VCON*  (,106/4)OFRF/2047 

NaRN(NRCTl)4l 

M«2 

IF  CNRC,LT,2)  GO  TO  11 


compress  DATA  OOMN  TO  AN  NCL  CELL  RADIAL. 


MaN»l 

00  1 K«2,m,2 
JlK/2 

00  1 I«l>3 

MVP(I»  JJ»(MVP(I,K*n+MVF(I,lO)/2,0 

continue 

n«nCL*4 

Ma2 


COMPUTE  DBZ 


DO  41  J»m,n 

PbmvP(3i JtEi 

IF  (P.LE.XCUT)  go  to  21 


USE  LINEAR  calibration. 


p«p«calm^calb 
GO  TO  31 


USE  NON  LINEAR  calibration. 


P«AC(nTAC(8)*P4AC(31*P**2*AC(4)#P**3 

w(J)«CKTP*20.*AL0G10  (SCON*  (FLOAT  (J«n*j5))*CELWTH(NTPtl)T,5 
IF  (w(J)  .LT.TLd))  w(Jl«0 


compute  V 


000004lj 
*000004lj 
0000041 
0000041 
0000041 
*0000041 
00000411 
00000411 
•0000042 
,0000041 
0000042, 
0000041 
0000042 
0000042 
0000042 
0000042 
•0000042 
0000042 
000004SQ 
0000043^ 
0000043 
0000043 
00000434 
0000043 
00000431 
0000043 
0000043 
00000434 
00000444 
0000044 
0000044 
0000044 
00000441 
00000441 
00000441 
0000044 
00000441 
00000441 
000004SC 
0000048 
0000048 
0000048 
00000481 
00000481 
00000481 
0000048 
00000481 
0000048 
00000441 
0000044 
0000044 
0000044 
0000044' 
0000044 
00000441 


70 


IF  (kJCJ)  .LE.TLd))  00  TO  41 
V(J)bIFIX(VC0N*FL0AT(MVF(1, 
8V(J)4lFlX(VC0N**2#Ft0AT(HVF(a, J*2))) 
IF  (VB(J) •EQ,«999,0R.NA.eQ.n  60  TO  41 
RB8C0N*{FLnAT(J.n.,5)*CEL«TH(NTF*l) 
V8(J)vCV(J)*VB(J))/R*1000. 

41  continue 

00  51  J«1,NCL 
31  VB(J)RVCJ) 

DO  61  JiliNCL 
61  VJ(J5«V(J) 

RETURN 

END 


00000467 

00000468 

00000469 

00000470 

00000471 

00000471 

00000473 

00000474 

0000047S 

00000476 

00000477 

00000478 

00000479 


o o r»  o o 


SUBROUTINE  prang 


VERSION  8,0  level  TS1119 
JHW  CDC6B00 

COMPUTES  RANGES  AND  PRINTS  THEM  OUT  FOR  BSCAN  MAPS, 


DIMENSION  R8AVE(S) 

C0MM0N/XNSUB/TL(0«) ,LT,T0M,DN,8TARTR,DELTR,RN(«) ,8C0N,CEL 
COMMON/ AD AT A /ID A V, I HOUR, I MIN, I8EC, NTP, N8F, NDD, NRC 

SCRAaSCON 

IF(NRC,ED,1)SCRA«8eON/2 

RMAX»8CRA*(RN(NRC*l)»,5)#CELWTHfNTP4n/lOOO, 

DaRMAX/8,0 

PSAVE(8)aRMAX 

JRT 

00  10  Iil,T 
R8AVE(J)«RSAVE(J«n«D 

10  continue 

CALL  PAGE 
WRITE(fe,99)R8AVE 

99  F0RMAT(IH0,31X,20HRANGE  SCALE  C KM  )/ 

XaX,2HAZ,«X,2HEL,lK.’3H0AV,lX,«HHHHM,lX,8HSS,8X,8FB,l,10H 

RETURN 

END 


wTM(3) 


OOOOOABO 
**•*00000481 
00000482 
00000483 
00000484 
****00000485 
00000488 
) 00000487 

00000488 
...•00000489 
00000490 
00000491 
00000492 
00000493 
00000494 
00000495 
00000498 
00000497 
00000498 
00000499 
00000500 
00000501 
00000502 
PRP)OOOOOS03 
00000504 
00000505 


o o r»  fi  o o o o ooooo 


SUBROUTINE  CONTOR  ( JM  Ax , KMAX,  IEM«X » I«T,'NIDi  NFC*  NZP«  NZH»  NUR*  NURU  K .'NOOOOOSOS 
iMZiNVl,NPA,lPvRNG,IPRNG,I0C,IDVC*IPBl»IPB8,IP83,IPCl*IPC2,IPC3,lPT0000050T 
2B,TB,1PBNT, IPCNT*T,IPTC.UP,TC# 1B»IC,HZ,VI, ICVNT, IBVNT,ni , TVB, IPTvBOOOOOSOS 
3*  lPVi»lPV2,lPv3,CI,  ATR,ZH,D8I,IPL0»TATR,KDD,NTT,C11,  IC1,NI»B  00000509 
R I IPBVNT* lACr* IACV*NUV*NVhAX*UV*VATR(NNE*BI)  00000510 


VERSION  2,0  level  TblllR 
JHW  AFGL  CDC6600 

fixed  contours*  perk  detection*  event  association. 


common  /ST0R2/  IMX 

COMMON  /InSUB/  TL(R) #LT,TDW,ON,STARTR*nELTR*RN(«) *8C0N*CElNTH(3) 
COMMON  /PARM/  PRlNTl  ,PRINT2*PRINT3,PRInT«*  ICODES  (3fc)  * A1  * B1  , A2, 82,' 
10NTRZ,C0NTRV,NFILE,‘nUMF,NREC,NUmR 
COMMON  /A102R/  MVPfS.lOZa) 

COMMON  /AZM/  AZMUTHCa«*0)  ,NA,ELEvAT,PRF;kEEP 
COMMON  /VALMAX/  ZMAX*VMAX,AC(«), CALM, CALB*XCUT,CK, INC 
COMMON  /AOATA/  inAy,lHniJR,lHlN,18EC*NTP,N8F,NDD,NRC 
COMMON  /MOREO/  INPRF,8CALE*L0V,LTV 
COMMON  /STORE/  AE* AA,B8,8L,CL,TV*T8V 

C0MM0N/w0RK/w(5i«),vC5ia5,vsf9ia),8Vf5ia),v8(5ia),vj(5ia) ,iii(5iai 
1 V8I (51«) ,H8 tsia) ,HV8JS1«5 ,NCL 
COMMON  /AZ2/8INA,Cn8A,0ELTAZ*I8CANF,NEL 

DATA  RPO/, 017453/ 

DATA  IPUP/3000/* IP0N/2nO0/ 

NTP»CELL  WIDTH  0,1,2  MEANING  ,5,1,042,2. 

NSFiSUBFRAmE  0. 1*2*3 

ndomfreq,  of  dump  pulses  ALT«0*ALL«1 

NRCRNO,  range  CELLS  0*1*2*3  MEANING  256,512*768,1024 

IDSLOTmO 

IF  (,NnT,PRINT3)  60  TO  1 
IF  (NA,E0,1)  call  prang 
IF(,nOT,CONTRZ)RET(IRN 

IF (NA,EO.n WRITE (2i  RN ( NRC+l) , CELWTH (NTP*!) , ELE V AT 
call  PRN3  f2,w) 


r00000511 

00000512 

00000513 

00000514 

>00000515 


LOGICAL  PRINTl, PRINTS, PRINT3,PR1NT4,C0ntRZ,C0NTRV  00000516 

integer  TCNTT) ,K00fNFCl ,IPvRN6(JMAX*IEMAX) ,IPRNG(JMAX,IEMAXJ ,IDCfI00000517 
IEMAX) , lOVCtlEMAXJ ,IPC1 ( JMAX,KMAx*lEMAXi , 1PC2 ( JM AX , KM AX , IEm AX) , IPC300000518 
2(JMAX,KMAX,IEMAX),!PTBfIEMAX),lPTC(IEMAX),IPBNT(KMAX,IEMAx),IPCNT(00000519 
3KMAX, IEMAX) ,IB(NPA, IEMAX,NFC) , IC (NPA, lEMAX* NFC) , TC (KMAX, IEMAX) , IPL00000520 
40(JMAX,KMAX) ,TB(KMAX, IEMAX) ,I8VNT(NFC) ;iCVNT(NFC)* ICl (NPA, IEMAX, NF00000521 
5C) , IPBI (JmAX,KMAX, IEMAX) , IPBP ( JMAX , KMAX , IEMAX) * IPB3 ( JMAX , KMAX , IEMA00000522 
6X) , lACT(NIO) ,IACV(Nin) ^IPVl C JMAX, KMAX, IEMAX) , IPv2 (JMAX , KM AX , IEM A X) 00000523 
7, IPV3( JMAX, KMAX, IEMAX) , TVS (KMAX, IEMAX)  , 'IPBVNT (KMAX* IEMAX) ,1PTVB(1E00000524 
SMAX)  00000525 

integer  W,V, VS*8V, VB, VJ,U1. VSI,H8,HVB  00000526 

INTEGER  TL*8TARTR,TV,T8V  00000527 

REAL  UP(NUP,NID) ,HZ(NHZ*NZH, IEMAX) , VI (NVi * NZH* IEMAX) * DI (IEMAX) ,08100000520 
1 (NID) ,CI(NPB, IEMAX,NFC) , ATR{IAT,NI0,NFC) ,ZH(nZP,NNE,NIO) ,TATR(NUmA00000529 
2X,NI0) ,CIl (NP8,IEMAX,NFC) , VATR(NVMAX,NID) *UV(NUV*NI0) , 00000530 

3 81  (NPB, IEMAX, NFC)  00000531 


•00000532 

00000533 

00000534 

C00000535 

00000536 

00000337 

00000538 

00000539 

00000540 

00000541 

00000542 

,00000543 

00000544 

00000545 

•00000546 

00000547 

00000548 

•00000549 

00000550 

00000551 

00000552 

00000553 

•00000554 

00000555 

00000556 

00000557 

00000558 

00000559 

00000560 


i 

I 


i 

j 

\ 


1 

continue 

00000561 

IF  (NA.Eo.ij  on  Tn  11 

00000561 

TEMPtAZMUTMINA.n 

0000056S 

4ZN0V(«AZPUTH(NA) 

00000564 

OCLTAZt(AZMUTH(NA)«TEMP)*RFD 

00000565 

TEMPiTEMP*RpO 

00000566 

80  TO  61 

00000567 

c 

0000056S 

c 

initialize. 

0000056B 

c 

00000570 

11 

TEMPaO.O 

00000571 

OELTAZrO.O 

00000572 

AZN0WR0,0 

OOOOOSTS 

00  21  KPl,NFC 

00000574 

21 

KOD(K)rO 

00000575 

SLi8IN(ELEVAT«RP0)/1000. 

00000576 

CLRCn8(ELEVAT*RP0)**2/AE/6,7SlE09 

00000577 

NCELPl 

0000057B 

NVCEL«1 

00000579 

ARpALOGIO(AA) 

000005S0 

RRr0.1/BB 

00000581 

ARiALOGIOCAA) 

00000588 

BRaO,l/BS 

00000563 

DO  31  KRlfNlD 

00000584 

DSI(K)«0, 

000005S5 

00  51  Jr1,NNE 

00000566 

00  31  L«1,NZP 

00000587 

11 

ZM(Lf JfKltO.O 

00000568 

00  «l  KP1,N10 

00000589 

DO  41  Jr1,IAT 

00000590 

00  41  L«1»NFC 

00000591 

41 

ATR(J»K,L)«0,0 

00000598 

DO  51  KRl.NFC 

0000059] 

DO  51  L«1*IEHAX 

00000594 

00  51  JR1,NPA 

00000595 

IB(J,L*K)«0 

00000596 

51 

XC(J>U.K)rQ 

00000597 

61 

CONTINUE 

00000598 

C 

00000599 

00  T1  K»1#IEHAX 

00000600 

DO  71  JR1,NZH 

00000601 

DO  71  Lil»NHZ 

00000608 

T1 

HZiLf J«K)rO.O 

00000603 

00  61  KRlflEMAX 

00000604 

00  SI  JRI.NVX 

00000605 

DO  SI  L«1»NZH 

00000606 

SI 

VltJ,L,K)«0,0 

00000607 

00  91  Kr1,IEMAX 

00000608 

DO  91  JBl.NFC 

00000609 

00  91  L»1*NPB 

00000410 

91 

CI(L*K,J)i0,0 

00000611 

DO  101  KRIiIEmAX 

00000618 

0I(K)»0,0 

00000613 

IOC(K)rO 

00000614 

T0vC(K)R0 

00000615 

74 


101 


111 


c 

c 

c 


131 

141 


C 

c 

c 

151 


c 

c 

c 

161 

171 


181 


191 


201 

211 


DO  101  jaitjHAX  00000616' 

IP9NG(JfK)«0  000006171 


IPVRN6(J,K)»0 
CONTINUE 
DO  111  KiliNFC 
ICVNT(K)«0 

continue 


00000618; 

00000619 

00000620 

00000621 

000006221 


00000623: 
00000628 
00000621 
00000626 
00000627 

FIND  events  00000628 

00000629 


IFV»0 

IFpO 

IFBaO 

IFVBPO 


00  281  l«2,NCl 
DO  231  Kil,NFC 
IF  (N(I),6T,TU(K))  60  TO  131 
00  TO  241 

IF  (Wd.n.LE.TLCKi)  GO  TO  Ul 
GO  TO  ISl 

ICVNT(K)PlCVNT(K)tl 

lEVENTilCVNT(K) 

IF  (K.EO.l)  lEOvIEVENT 

ICCl»IEVENT.K)«I.l 

ICC3*IEVENT,K)«IE0 

TAULT  ATTRIBUTES# 

R»aC0N*(FL0AT(I«l).,5)*CELWTHCNTF4l) 

lEVENTilCVNTtKJ 

CIfl,IEVENT,K)«CI(l.IEVENTftK)TR 

CI(2,IEVENT,K)«CIt2,IEVENT,K)+B*W(I) 

CIC3,IEVEnT,K)»CI(3,IEVENT,K)4-R*R»W(I) 

IF  (K.NE#n  GO  TO  231 

PEAK  DETECTION.  LOCATE  AND  COUNT  PEAKS, 

IF  (W(n.w(I«n)  ITl.lBl.'lSl 

IPBiI-1 

60  TO  181 

IF  (IPB.EQ.O)  GO  TO  181 
IP«IPT1 

IPRNG(!P,IEVENT)*(ItIPB)/2 

IPBpO 

continue 

IF  (VSd)  ,E0,»999)  GO  TO  191 

IF  (V8fl»n  ,E0,»999J  GO  TO  201 

IF  (IAB8(V8(I))«IAB8(V8(I«n))  l9l,2ir;201 

IF  (IPVB.EQ.O)  GO  TO  211 

IPWiIPVTl 

iPWRNGdPV,  IEVENT)B(lTlPVBi/2 

IPvBtO 

60  TO  211 

1PVBPI«1 

CONTINUE 


00000630 

00000631 

00000611 

00000633 

00000638 

00000631 

00000636 

00000637 

00000638 

00000639 

00000640 

00000681 

00000681 

00000683 

00000688 

00000681 

00000686 

00000687 

00000688 

00000689 

00000618 

00000691 

00000692 

00000693 

00000698 

00000691 

00000696 

00000697 

00000698 

00000699 

00000668 

00000661 

00000668 

00000663 

00000668 

00000669 

00000666 

00000667 

00000668 

00000669 

00000670 


75 


r 


221 

2S1 

2«1 


C 

C 

c 


251 


261 


2T1 

281 


» 

I 


( 

I 


291 


801 

311 

321 


IH«IFIV(I?*SL9R*R*CL)«1 
IR  (IH.LEtO.OR.IH.OT.NlH)  60  TO  221 
If l»IfVENT 

HZtl»XH,XEl)fM2(l,lH,lEn*WCI)*R 

HZ(2f  XH,IEl)«HZ(2,XH.XCn<»R 

XR  (W(1).6T.TV.0R.8V(X5,6T.TSVJ  60  TO  221 

XR  (w(X).LT,TL(n,0R.VCX).E0,»999)  60  TO  221 

VX(lilH,ICl)«VX(liXH.XEn«V(I) 

VX(2iXH,XEn«VX(2,IH,XEnRV(I)*V(X) 

VX(3|lH»ZCl)«Vl(S|XH,XEn9l.O 

continue 

XR  (NCLtNE.l)  60  TO  231 
RAZN«10,**rBR*W(Z)«AR} 

OX(lEl)»nX(XElWRAlN«R 
CONTINUE 
60  TO  281 
DO  271  KLiKfNRC 

XR  (W(X*1),LE.TL(KL))  60  TO  281 

iewent»xcvnt(kl) 

XCr2iIEVCNT,KL)«X»l 

KEEP  COUNT  OR  PEAKS  WITH  EVENT, 

XR  (KL.NB.l)  60  TO  271 
XR  (IP8.EQ.0)  60  TO  2Sl 
XP»XP*1 

XPRN6(XP,XEVENT)*(!*XPB)/2 

XPB«0 

XDCCXEWENTJiXP 

XPaO 

XR  (XPVB.EQ.O)  60  TO  281 
XPV»IPV81 

XPVRn6(XPV,XEVENT)b(X^XPV8)/2 

XPV8»0 

1DVC(XEVENT)«XPV 

Xpvao 

CONTINUE 

continue 

XR  (NA.NE.n  GO  TO  321 
DO  311  KaifNRC 
DO  311  KEVENTk1,XEHA¥ 

DO  291  Xk1,NPB 

CX1(X,KEVENT,K)«CI(X,KEVENT,K) 

DO  301  XalfNPA 

XC1(X#KIVENT,K1bXC(I,KEVENT,K) 

continue 

C0ST»C08(TEMP) 

8XNT«8IN(TEmP) 

CO8A«Cn8(AZNOW*RP0i 

C09A2«C08a*C08a 

SXNAaSXNlAZNONKRPOi 

8XNA2b8XNA*8XNA 

8NACNAa8ZNA*CQ8A 

XMyaRN(NRC9l)/2«3 

IRfNRC.Eq.nZMXBNRC»3 


00000871 

00000872 

0000087S 

00000878 

00000871 

00000878 

00000877 

00000878 

00000879 

00000880 

00000881 

00000881 

000008BS 

00000884 

00000885 

00000888 

00000887 

00000868 

00000889 

00000890 

00000891 

00000892 

00000891 

00000894 

00000895 

00000898 

00000897 

00000898 

00000899 

00000700 

00000701 

00000702 

00000703 

00000704 

00000705 

00000708 

00000707 

00000708 

00000709 

00000710 

00000711 

00000712 

00000715 

00000714 

00000715 

00000718 

00000717 

00000718 

00000719 

00000720 

00000721 

00000722 

00000725 

00000724 

00000725 


76 


j t 


PLOT  FIXED  CONTOURS. 


DO  *11  KSI.NFC 

IFU»IPUF*K 

IPO»IPON*K 

IDO«KDD(K) 

KEVENTll 

XEVCNTtl 

IF  (IB(2»IEVEN 

IF  (IBClflEVEN 

IF  (IB(2,IEVEN 


■ 1 

• 1 

(2»IEVENT,K)^E0.0.AND.IC(2iKEVENTfK) .EQ.O) 
(l»IEVENTfK),ST.IC(2/rKEVENT,K))  CD  TO  ATI 
(2,IEVENT,K),LT,IC(1,KEVENT,K))  go  to  ATI 


60  TO  *01 


ASSOCIATED 
LEFT  SIDE  PEN  UP. 

IIOAlB(NPA,IEyENT,K) 

IF  (IBCANF.EO.O)  IC(NPA,KEvENT.K)aIID 

XbFL0AT(I8(1,IEVENT,K))^.B 

Rb8C0N*X*CELWTH(NTPb1)/(3.8A*10E03) 

Xb8CALE*(R*8INT*A.O) 

ya8CALE*(R*C08TB4.0) 

IF  (PRINTa)  CALL  PLOT  (X,Y,3) 

WR1TE(2)X,V,IPU 

LEFT  SIDE  PEN  OONN. 

X»FL0AT(IC(1,KEVENT»K)5».'5 

Rf8C0NBXBCELHTH(NTPAn/CJt8A*10E03) 

XaSCAIE*(R*8XNA««.0) 

YbSCALE*(R*C08a«A.O) 

IF  (PRINT«)  call  plot  (X,Y,21 
WRITE(2)X,Y,IPD 

ATR(l,lID,K)»ATR(ltXID,X)BDELTAZ*CI(l,KEVENT;K) 

ATRCIf I1D.K)BATR(24IID(K)B0ELTAZ*CI(2.KEVENT«K) 

ATR(3f  IID,K)iiATR(3.IIO,K)bsXNA*OELTAZ*CI  (SiKEVENT.K) 
ATR(A,lID,K)iATR(A,IX0,K)»C08A*DELTAZACX t3,KEVENT,K) 
IllfIC(S,KEVENT,K) 

IIDllXCCNPA.IEl,!) 

IF  (ATR(IAT,IID,K):EQ.0.)  ATR(IAT,IID,K)bI101 
IF  (IC(1»KCvENT»K).E0.1.0R.XC(2,KEVENT.K).CQ.IMX)  atrcxat.iid, 
1ABS(ATR(IAT,IID>K)) 

IF  (X.NE.n  80  TO  371 
I2fNNE 

DO  3*1  IHbi.NZH 

IF  (HZ(2iIH,KEVENT)  .LE.'Ol  60  TO  3*1 
IF  (VI(3.IH.KEVENT) .LE.O)  60  TO  SSI 
ZH(1,I2iII0)bZH(1.I2«IID)pvX(1iIH,KEvENT) 
ZH(2,I2,X!D)nZH(2.I2*XI0)«VI(2.IH.KEVENT) 

ZH(3f Z2,IID)aZH(3,I2.II0)A8INA*vI(1»IH^KCVENT) 
ZH(A,It,IID)aZH(A.I2.IID)BC08AAVl(lf IH.KEVENT) 
ZHfS,l2,Xl0)BZH(5,I2,II0)P8XNA2BVlC3,IH,KEVENT) 
ZH(6»X2,IID)"ZH(*,12.1ID)PC0SA2*VI(3,Ih,KEVEnT) 

ZH(7>l2f XI0)«ZH(7.I2.IX0)A8NACNA*V1(3,IH»KEVENT) 

TH(8,12f XI0)AZH(6,I8.XID)«8INA«vX(3»IH.KEVCNT) 
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OOOOOTSO 
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0000074* 

00000747 

000007*8 
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351 

361 

371 

381 

C 

C 

C 


C 

C 

C 


C 

C 

c 


401 

411 

C 

421 


ZH(9,I2,II0)«2H(9,I2 

ZHflO,I2,nD)«ZH(10, 

ZWfll»I2,110)«ZHCll, 

ZH(12,X2,II0)«ZH(12, 

COMTINuC 

08I(II0)b08X(XID)«0X 
IF  (ISCANF.CO.I.AND. 
XF  (XB{1,IEVENT41,K) 
XF  (XB(1, IEvENT^IiK) 


,XXO)^COSA*vX(S,XH,KEVENT) 

la.XIDlA'VltSiIHiKEvENT) 

X2,XI0)*HZ(1,XN,KEvENT)*0ELTAZ 

I2,XIDi*HZ(2iXH,KEVENT)*0CLTAZ 

(kevent)*oeltaz 

XlR,NE.XC(NpA,KEVENT,Kn  GO  TO 
,GT,ICf2,KEvENT,K))  GO  TO  44l 
•EQ.O)  GO  TO  441 


1 


571 


DRAW  DOWN  TO  PRESENT  AZMUTH. 

XpFL0AT(X8(2,XEvENT,K))«,5 

RB8C0N*X*C€LWTH(NTP4n/(S,R4*lOE03) 

Xb8CALE*{R*8XNT*4,0) 

Y»8CALE*CRaC08T44,0) 

IF  XFRXNTa)  call  PLOT  fX»Y,3) 

WRITE(2)X,V,XPU 

Xp8CALE*(R*8XNA«4,0) 

VaSCALE*(R*C08A«4.0) 

IF  (PRXNT4)  CALL  PLOT  (X,Y.2) 
wRITE(2)X, Y, IPO 
WRITE(2)X,V,IPU 

DRAW  OVER  TO  lEVENT*! 

X«FLnAT(lB(l,XEVENT*l,X)l«,5 

R«8C0N*X*CELWTH(NTP*n/(3,«4*10E03) 

X«8CALE*(R«SXNAP4,0) 

Yp8CALE*(R*C08A«4,0) 

IF  (PRXNT4)  CALL  PLOT  fX.Y/E) 

WRITC(2JX,Y,XP0 

WRrTE(2)X,Y,IPU 

DRAW  UP  TO  PREVIOUS  AZMUTH, 

Xa8CALE*(R*SXNTP4.n) 

Yp8CALE*CR*C0ST*4,0) 

XF  (PRXNT4)  CALL  PLOT  tX.Y.’B) 

WRITEC?)X,Y,XP0 

IEVENTiIEvEnT+I 

XF  flEvENT.GT.IEMAX)  GO  TO  601 
KXn*lB(NPA, lEVENTjK) 

IF  {ATR(lAT,KXD,K),E0,0,0,0R,ATR(IAT,XI0,K) .EO.O.O) 

IA7T«XAT*1 

DO  411  Jtl.IATT 

ATR(J»I10,K)pATR(J,1X0;k)*ATR(J,KI0,K) 

XF  (ATR(XaT,KXD,K)  .’LT,0,0,0R,ATR(XAT,XID,K)  .LT.O.O) 
1««A88CATR(XAT,IX0,K)) 

A 0 will  flag  useless  attr’b, 
ATRtIAT»KI0,K)»0,0 
IOBLOTpXID 
I2«NNC 

ZH(J,X2,II0)PZH(J»X2.II0)tZH(J,l2,KXD) 

ZHfNZP.IEfRIOlwO.O 


78 


GO  TO  381 
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00000619 


F 


osi(iiD)iosi(iio}«osr(Kro) 

Gn  TO  1S1 

a«l  IF  {IC(l»t<EVENT»l,«)  ,CT,I0(2,'IEvEnT,K))  GO  TO  «5i 
IF  CIC(l,KeveNT*l,F) ,E0,0)  60  TO  «5l 

c 

C DFAW  LXWt  C0N»JECTInG  IC(n)  TO  IC(N*l), 

C 

*«FL0AT(IC(2»FEVENT,K))-,5 

F«scoN*x*ceLWTH(MTF*n/(3.sa*ioEos) 

X«SCALE*(R*Sl^A^a.O) 

y«8CALE*(R*C0SA«4,0) 

IF  (PRINT#)  CALL  PLOT  fx;T,3) 

WRITC(2)X,Y,1PU 

XiFL0AT(lC(l»KCVENT4liK))-,5 

R«SC0N*X*CELWTM(NTP*n/(3.ia*10E03) 

XbSCALE«(R*SINA««,0) 

V9SCALe*(R*C08A«4,0) 

IF  (PRINT#)  CALL  PLOT  (X,r,2) 

wRITE(2)X,Y,IP0 

KEVENT#KEVCNT*l 

IC(NPA,KCVCnT,K)«IIO 

IF  (KEVENT.LT.IEMAX)  60  TO  3#1 

60  TO  801 

C 

C RIGHT  SIDE, 

C 

#51  X*FL0AT(IB(2,IEVENT,K))«.5 

R«8C0N*X*CELHTH(NTP*1)/(3.8#*10E03) 

XaSCALE*(R*8lNTt#,0) 

Ys8CALE*(R*C0ST^#,0) 

IF  (PRINT#)  CALL  PLOT  (X.Y.S) 

WRITE(2)X,Y,IPII 

X«FL0AT(IC(2,KEVENT,K)).,5 

R«8C0N*X*CELNTH(NTP*1)/(3,8#*10E03) 

XI8CALE*(R*8INA##.0) 

YaSCALE*(R*C08A«#,n) 

IF  (PRINT#)  CALL  PLOT  fX,Y,2) 

WRITE(2)X,Y,IP0 

lEVENTilEvENT^l 

IF  (IEvENT.GT.IEHAX)  go  TO  601 

KEVENT«KEVENT*1 

60  TO  331 

#71  IF  (I8(1,IEvENT,K),E0,0)  60  TO  521 
IF  (IC(1,KEVENT,R)  ,EQ.O)  GO  TO  #81 

C 

C UNA880CIATC0. 

C angle  LINE  ON  IB 

C 

IF  (IC(2,KEvENT,K),’lT,IB(1,IEVENT,K))  so  to  511 

#81  X«FL0AT(IB(1,IEVENT,K))-,3 

R«8C0N*X<»CELNTH(ntp*1)/(3.8#*10E03) 

X«SCALC*(R*8lNTt4,0) 

YR8CALE* (R*C08T##,0) 

IF  (PRINT#)  CALL  PLOT  (X.’Y.'S) 

WRITE(2)X,Y,IPU 
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00000843 
0000084# 
00000845 
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00000847 
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00000846 
00000850 
00000851 
00000852 
00000853 
00000854 
00000855 
00000856 
00000857 
00000858  1 

00000856  I 
00000860  : 
00000861  > 
00000862  ; 
00000863 
00000864  j 
00000865  ! 
00000866 
00000867 
00000868  ! 
00000866 
00000870 
00000871  ' 
00000872 
00000873 
00000874 
00000875 
00000876 
00000877 
00000878 
00000876 
00000880 
00000881 
00000882 
00000883 
00000884  I 
00000885 
00000886  1 
00000887  i 

00000888  I 
00000886 
00000860 
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XSvICALC* 

VaSCALr^C^^CO^A^tt.O) 

If  (PRI^T«)  CALL  PLOT  tXfV.21 

wRiTec2)*,Y,iPn 

WRITE(*) *, V, IPU 

X»PLnAT(IB(2,IEvENT,K)).,5 

R«8C0N*X*CELWTH(NTPtl)/t1.8«*l0C03) 

»«SCALE*CR*8IMT^a,0) 

ViSCALE*(R*COST«a,0) 

IF  (PRINTa)  CALL  PLOT  f*,V,21 
WRITE(?.)X,  Y,  IPO 
1I0»18(NPA,IEVENT,K) 


lEl»lCt3,KEVENT,K) 

IID1«ZC(NPA,IE1,1} 

IF  (ATR(IAT,IIO,K),EO,0.)  ATR ( I AT, I ID , K ) ■ 1 1 D I 
IF  (IC(1,KEvENT,K) ,EQ,1,0R.IC(2.KEvEnT,K) .EO.IM: 
1A88(ATR(IAT,IID,K) J 
lEVENT-IlVENT^l 

IF  (IEvENT.CT.IEMAX)  go  to  601 

If  (IC(i,KEvENT,K),LE,IB(2,IEVENT,K))  CO  TO  33l 
IF  (IC(2,KEvENT,K)  ,NE,0)  GO  TO  501 
GO  TO  331 

)l  IF  (IB(1,IEVENT,K),EO,0)  go  to  5Pl 

.1  IF  tIC(l,KEvENT,K),GT,IB(2,IEVENT,K))  CO  TO  Hi 

UNA8S0CIATE0 
STRAIGHT  line  ON  IC, 

•A  IP  (IC(l,KEVENT,K)  ,EQ.O)  GO  TO  562 
IF  (ISCANF.EQ.n  GO  TO  581 
IFtID8LOT,EO.O)GO  TO  523 
ICfNPA,KEVENT,K)»IOSLOT 
GO  TO  52tt 
523  I00»ID0+l 

ICfNPA,KEVENT,K)«IOD 
IF  (NA.EQ.n  ICltNPA,FEVENT,K)«lOD 
52a  X»FLOAt{IC(1,KEvENT,K)5»,5 

R«8CnN*X*CELWTH(NTP^l)/(S.8a*lOE03) 

*«8CALE*(R*SlNA»a,0) 

YaSCALE* (R*C08A«tt,0) 

IF  (PRINTaj  CALL  PLOT  (X,Y,3) 

WRITE(2)X,Y,1PU 

XaFL0AT(lC(2,KEVENT,K))*,5 

RaSC0N*X*CELWTM(NTP^n/(3.8a*10Ee3) 

XaSCALE*(R*8lNA44,0) 

Ya8CALE*(R*CnSA«a.O) 

IF  (PRINTa)  CALL  PLOT  (X.Y.’E) 

WRITE(2)X,Y,IP0 
IF  (na.EQ.I)  go  to  531 
IF(ID8L0T,E0,0)G0  to  52T 
lOTEMPalDO 


ATR(IAT,TID,K) 
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0000083 

0000083 

0000083 

0000083 

0000083 

0000083 

0000088 

0000084 

0000084 

0000084 

0000084 

0000084 


80 


IDOpIOSLOT 

00000998 

527 

ATB(l,IDR,K)PnELTA7*CI(l,KEVENT,K) 

00000997 

ATR(2, I00,K)«DELTAZ*CI (2,KEVENT,K) 

00000998 

ATR(3,ID0,K)bSINA<»DELTAZ*C1  (S,KEVENT,K) 

00000999 

ATRCa,IDD,K)«C08A*0ELTAZ*CI (3,KEVENT,K) 

00000950 

531 

1EIPIC(3,KEVENT,K) 

00000951 

ATR(IAT,IDD,K)pIC(npa,IE1,1) 

00000952 

IP  (NA.EQ.n  ATR(IAT,IOO,KJp«ABS(ATR(IAT,IDO,K)) 

00000953 

IP  (IC(1,KEVENT,K)  .EO.l.OR.lCta.KEVENT.'K)  .EO.IMX) 

ATk Cl AT, IDD,K) ••00000959 

IAB8(ATR(IAT,I0D,K)) 

00000955 

IP  (P.NE,!)  60  TO  581 

00000958 

I2PNNE 

00000957 

no  551  IH«l,NZH 

00000958 

IP  (HZ(2,1M,KEVENT) ,LE,0,)  60  TO  551 

00000959  1 

IP  (Vl(S,IH,KEVENT),Ea.O.)  60  TO  591 

00000980  . 

ZHfl,12,IDD)PVl(l,IH,KEVENT) 

00000981  ‘ 

ZH(2,12,100)pvI(2,IH,KEvENT) 

00000982  ; 

ZH(3,12,  IDO)  PSINAovI  (!,' IH,KE  VENT) 

00000983  ! 

ZH(a,I2,IDD)PC0SA*Vl(l.IH,KEVENT) 

00000989  I 

ZK(5,I2, ID0)P8INA8*VI(1,IH,KEVENT) 

00000985  ! 

ZH(8,I2,I00)aC0SA2*vl(3,XH,KEVENT) 

00000988 

ZHtT, I2,lDD)P8NACNA*Vlf3,lH,KEVENT) 

00000967 

ZH(8,I2,I0D)PSINA*vl(3,IH,KEVENT) 

00000986 

ZHf9,I2,ID0)PC08A*Vl(1,lH,KEVENT) 

00000969 

ZH(10,I2,ID0)pv1(3,IH,kEVEnT) 

00000970 

591 

IP  (NA.EO.n  GO  TO  581 

00000971 

ZH(11,I2»I0D)pHZ(1;IH,KEvEnT)*0ELTAZ 

00000972 

ZHn2»I2,I00)pMZ(2,IH,Kev£NT)»DELTAZ 

00000973 

551 

CONTINUE 

00000979 

08I(IDD)pDI(KEVENT)*DELTAZ 

00000975 

581 

IPdOSLOT.NE.O)  IOOpIOTEMP 

00000976 

IPSLOTpO 

00000977 

582 

KEVENTpKEvENT^I 

00000978 

It  (KEVENT.GT.IEMAX)  go  to  801 

00000979 

GO  TO  331 

00000980 

571 

KI0*IC(NPA,KEVENT,K) 

00000961  I 

IP  CICd  •KEVENT.K)  .EO.1.0R,IC(2,KEVENT,'k)  .EO.IMX) 

GO  TO  901  00000982  1 

ATP(IAT,Kin,K)pA8Sf ATR(IAT,KID,K)) 

00000983 

GO  TO  ttOl 

00000969  ; 

581 

IIOpIC (NPA,KEVENT,K) 

00000985  f 

591 

IP  CICfl,KevENT,K) ,0T,1 ,AND,IC(2,KEVENT,K) ,LT,ImX) 

ATR(IAt,II0,K)«00000988  ' 

1A88(ATR(IAT,II0,K)) 

00000987 

IP  (ATR(IaT,IIO,K),NE.O.)  go  to  391 

00000988 

GO  TO  591 

00000989 

801 

Knn(K)*IDO 

00000990 

811 

CONTINUE 

00000991 

782 

IP C,NOT,CONTRV)Gn  TO  800 

00000992 

CALL  PEAKO  ?W,LDV,TLfn.IPRNG,IOC,l,TATR,IPBl,!PB2 

, IPBS,IPTB,TB,IP00000993 

IPCNT,T,1PC1 ,IPC2, IPC3;1PTC,UP,TC,NTT, 1EMAX,KMAX* JMAX,  NCE0000099a 

2L,''I0,IB,  IC,  XCVNT,  I8VNT,^iPA,NPC,  IPLO,Nl)MAX,NLlP»  IACT»H8)  00000995 

C*UU  PEA<<0(V9,LTV,O,IPV9n0;I0VC,0,  VATR.'IPVI,  KPV?,  IPV3,  00000998 

*I»Tv9,TVB, IPBVNT, IPCMT.T, IPCI , IPC2 , IPCI , IPTC , UV , TC , NTT , IEm AX , KM AX , 00000997 
1J-AI,mvCEL,nI0,1B,IC, ICVNT.IBVNTjNPA.NFC, IPL0»NVMAX,NUV, IACV,HV8)  00000998 

00000999 

STORE  PRESENT  PARAMETERS  IN  PREVIOUS  PARAMETERS,  00001000 


81 


n 


Name»41 

00001431 

WRITE (6,9909) NAME. TNDX 

00001435 

365 

continue 

00001440 

GO  TO  931 

00001441 

C 

00001442 

C 

combine  NPCEL  ANO  LPCEL,  PEAK  VALUES  EQUAL 

00001443 

c 

00001444 

c 

00001445 

c 

COMBINE  WITH  B RADIAL  CELLS 

0000144* 

c 

00001447, 

421 

IF(mPK,L7,0)G0  TD  422 

00001448 

lpcel«mpk 

0000144* 

IF(TATR(I0X, LPCEL) .EO.NA, AND, NPK,EQ,D)fiO  TO  485 

00001450 

IN0X»T AT»(1, LPCEL )*TC(KC, IE) *1 

00001451 

IFCINDX.LT.OGO  to  481 

00001452 

c 

00001453 

C 

COMBINE  WITH  B • RADIAL,  C-LEVEL  LOwER 

00001454 

c 

00001455 

IPC3(1PE,KC,IE)»LPCEL 

00001458 

iNlINOXALM 

00001457 

512 

ISTBlPCl (IPE.KC, lEiwl 

00001458 

ISP«IPC2(IPE,KC,IE) 

00001459 

NPCEL«LPCEL 

00001460 

DO  511  I«I8T,I8P 

00001461 

R«SC0N*(FLnAT(I«l)*,5)wCELWTH(NTPTl) 

00001462 

TATR(2t1N, NPCEL )«TATR(2tIN,NPCEL)+0AZ*R 

00001463 

TATR(3tIN,NPCEL)»TATR(1tTN,NPCEl)*DAZ»R*U(I) 

00001464 

TATR(4*IN,NPCEL)«TATR(4TlN,NPCEL)T0AZ(»8AZ*R*RWU(t) 

00001465 

TATK(5wIN,NPCEL)»TATR(5TlN,NPCEL)TDAZ*CAZ«ifRwR*U(I) 

00001468 

IF  (ITV.NE.i)  go  to  531 

00001467 

IF  (V(I)  ,Ea,»999,0R,V(I»l) ,EQ,-999)  GO  TO  521 

00001468 

TATH(6TIN,  NPCEL)  ■TATR(6tIN,  NPCEL)tDAZ*R*(V(I)**V  (1-1)) 

00001469 

521 


5Jl 


IF  (VSfl)  ,E0,*999)  Gn  TO  5S1  OOOOUTO 
TATB(7»IN,nfCEL)«T»TR(T  + IM,NI»CEu)*DAZ*B«>V8(I)  0000  t«Tl 
TATRt8*lN,NPCEL)«AMAXl (T ATP («+lN , NPCEL j . FLO AT ( I ABS ( V* Cl) ) ) ) 00001472 
continue  00001475 
TATR  (ir»X*IN,NPCEL)  ■SIGN  (FLOAT  (NA  ),  TATS  CIDX  + IN,NPCEL)  ) 00001474 
IF  (I8T,E0.2.0R,1SP.EO,IMX)TATR(10X*IN,NPCEL)»SIGN(TATR(IDx  + IN,wPCE00001475 


C 

C 

c 

c 

c 

C 


461 


XL). •1.0) 

NAHEaSl 

write (4,9909)NAHE, in 

WRITECS,  IOTDNPCEL.  (TATRCKZ.NPCEL)  .KZ«1,NUWP) 

GO  TO  422 

COMBINE  with  a. radial.  c«level  higher 

IF  FIRST  COMBINE.  AREA.O,  IF  SECOND  OR  HIGHER, 
TEST  AREA  TO  ESTABLISH  NEW  NUMBERS 

INDX»»IN0X 

INO^NUWP 

lNSi2 

TATR(l,LPCEL)«TC(KC.IE)*l 
TATR(NUMP,LPCEL)«IC(NPA.IE. 1) 

IF(INDX,GE,LDB)GD  to  482 


AREA!*! . 


000014TB 

00001477 

0000147a 

0000147a 

00001480 

000014S1 

000014S2 

00001485 

00001484 

00001485 

0000146* 

00001487 

00001488 

00001485 

00001490 

00001491 

00001491 
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r 


800  DO  790  KiliMFC 

DO  790  IEVENT»l,ieMAX 
DO  790  N«1.NP6 

790  8I{N,IEVENT,K)«CICN,IEVENT,K) 

DO  801  KvifNFC 
IBVNT(K)>1CVNT(K) 

00  801  lEvENTul.IEMAX 
00  801  N»l,NPA 

IBrN,ZEVENr,K)«ZC(N,ievENT.K) 

801  ZC(N»ZEVEnTiK)«0 

ZF  (ZSCAMF.EO.l)  80  TO  871 
RETURN 
C 
C 

ENTRY  C0NTR2 

ZF  (18CANF.GT.0)  60  TO  831 
DO  821  KRliNFC 
lEaZBVNTCKl 
DO  821  ZilfZE 
ZDD«ZB(nPA,ZE,K) 

821  ATR(ZAT,Z0D,F)»-ATR(ZAT, Z00,K) 

C PLOT  FZNAL  RA0ZAL8, 

X*4,0 

v«a,o 

ZF(PRZNTa)CALL  PL0T(X,Y,3) 

wRZTE(2)X,Y,ZPU 

SCRA«8C0N 

IF(NRC,EQ.3)8CRA«8C0N/2 

R«SCRA*(RNtNRC*n»,5)*CELWTH(NTP*l)/(3,’84«H0E0S) 
X»SCALE*(R*8ZN(AZMlJTH(n*RPD)t4,0) 
Y«SCALE*(R*COS(AZMlJTH(n  *990)  44,0) 

ZF(PRZNT4)CALL  PL0T(X,Y,2) 

WRZTE(2)X,Y,ZP0 

X»8CALE*(R*8ZN(AZHUTH(NA)*RPD)t4.0) 
Y»8CALE*(R*C08(AZmuTH(NA)*RP0)^4,0) 
ZF(PRZNT4)CALL  PL0T(X,Y,3) 

WRlTE(2)X, Y, ZPU 

XB4.0 

Ya4,0 

ZF(PR1NT4)CALL  PL0T(X,Y,2) 

WRITE(2)X,Y, ZPO 
GO  TO  871 

831  DO  861  KHfNFC 
ZE«ZCVNT(K) 

DO  861  ZaiiZE 
DO  841  L>liNPA 
841  ZC(L,l,K)aZCl(L,Z,K) 

DO  891  LRlfNPB 
891  CZ(L»liK)aCZl(L»Z,K) 

861  CONTINUE 

TEHP«AZMUTH(NA)«RPO 

OELTAZ»fAZMUTH(n*TCMP)*RPD 

AZN0W»AZMUTH(1) 

80  TO  321 


00001001 
00001002 
00001003 
00001004 
00001009 
00001006 
00001007 
00001008 
00001009 
00001010 
0000101 1 
00001012 
00001013 
00001014 
00001019 
00001016 
00001017 
00001018 
00001019 
00001020 
00001021 
00001022 
00001023 
00001024 
00001029 
00001026 
00001027 
00001028 
00001029 
00001030 
00001031 
00001032 
00001033 
00001034 
00001039 
00001036 
00001037 
00001038 
00001039 
00001040 
00001041 
00001042 
00001043 
00001044 
00001049 
00001046 
00001047 
00001048 
00001049 
00001090 
00001091 
00001092 
00001093 
00001094 
00001099 


I 
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n o c»  r»  r»  o « ^4  usfeor»o»<i  r»r»oo 


pun*  c event  loop 

DO  951  lE«l,Xe*^ 

1PL»1 

IPtIDCdE) 

IF  (IP.EQ.O)  GO  TO  951 

Jll«0 

JESaO 

FIND  6 EVENTS  ASSOCIATED  WITH  C EVENTS, 
JEM  IS  NO,  OF  EVENTS  IN  PREVIOUS  RADIAL, 

JEMalBVNTCl) 

IFtJEM  ,EO,  0)  SO  TO  ai 
DO  31  JC«1,JEM 

IF  (tC(«.lE,n.NE.IB(9;jE,t))  GO  TO  31 
JE2PJE 

IF  (JEl.EQ.O)  JElaJE 

continue 

find  THRESHOLDS  FOR  IE  EVENT 

DO  51  J"l,JMX08 
T(  J)aO 

DO  Tl  LPlPL,IP 

IRiIPCRN6(L,IEl 

00  71  Kil,LOB 

ITalAB8CU(IR)5*TM«K4l 

IFtlT,GE.l,AND,IT,LE.JMXOB)TfIT)ai 

continue 

IPTal 

DO  91  L«l,JMXDB 

IF  (T(L))  91,91,81 

TC(lPT,IE)iL*TM»l 

lPTPlPT+1 

CONTINUE 

IPT»IPT«1 

IPTCfIE)iIPT 


LOOP 


RANGE 


IE  EVENT  TO  FIND  CONTOUR 


IB6N«IC(1,1E,1)«1 
IN0aIC(2,lE,l)^l 
DO  161  IalBGN,lNO 

IF  (I,NE,IPCRNG(IPK,IEn  SO  TO  lOl 
IPKBIPK#! 


LOOP 


threshold 


131  K«l,lPT 

(u(i) ,eo,«999)  go  to  lai 
(lA0S(U(n),GT,TC(K,IE))  Gl 
TO  lai 

(U(I«1) .EG.^RGR)  GO  TO  121 


00001218 

000012181 

0000122a 

000012211 

00001228 

0000122S 

00001224 

00001225 

0000122 

00001227 

0000122 

0000122 

0000123 

00001231 

00001212 

00001233 

0000123 

00001235 

0000123* 

00001237 

00001238 

0000123* 

0000124 

00001241 

00001242 

00001243 

00001244 

00001245 

00001248 

00001247 

00001248 

0000124*! 

00001250 

00001251 

00001252 

00001253 

0000125*! 

00001255 

00001258i 

00001257 

00001258 

0000125* 

00001260 

00001261 

00001262 

00001263 

0000126* 

00001265 

00001268^ 

00001267 

00001268 

0000126* 

00001270 

00001271 

00001272 
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oor»  ooooo  r>  nonoo 


t 


11 


21 

C 


SUBROUTINE  PEARO  (U, LOB,  TM,' IPCRNG,  IDC,  ITY , TATR,  IPBl , IPB2,  XPB5,  IPTBOOOOl  U3 
1,TB,IP8NT,IpCnT,T,IPC1,IPC2,IPC1,IPTC,iiP,TC,  JMXOB,IEMAX,KmAX,  jMA>r,0000116a 
2NCELLfNlO,IB,lC,ICVNT,IRVNT,NPA,NFC,IPLO>NUHAX,NUP,IACT,HB)  00001165 


VERSION  !• 


0 LEVEL  T70112 
JHM  AP6L  6600 

DETERMINES  PEAK  VALUES  AND  THEIR  ATTRIBUTES, 


•00001166 


00001167 

00001166 

00001169 


•00001170 


integer  IPCRNGtJMAX, IEMAX) , lOC(lEMAX) , IPBl ( JM AX , KMAX , lEMAx) , IPB2 f JOOOOl 1 7 1 
1MAX,KMAX,IEMAX1 ,IP85(JMAX,KMAX,IEMAX1 ,1PTB(IEMAX) ,TB(KHAX,IEMAX) ,100001172 
2P8NT(KMAX, lEHAX) ,T(JMX0B) ,IPC1 (JMAX,KHAX, IEMAX) ,1PC2(JMAX,KMAX*IEM00001173 
3AX) ,IPC3(JHAX,KMAX,IEMAX) ,IPTC(IFMAX),TC(KMAX»IEMAX) ,IPL0(JMAX,KMA0000117« 
RX),IB(NPA,IEHAX,NEC),IC(NPA,IEMAX,NFC);iCVNTfNFC),IBVNT(NPC),IPCNT00001175 
5 (KM AX, lEH AX), I ACT  (N 10) ,H8(1) ,U(l)  00001176 

integer  W, V, VS,SV,VB, VJ,UI. VSI,H1 ,H2  00001177 

integer  TV,T8V,TM,TL,8TARTR  00001176 

real  TATR(NUMAX,NI0) ,UP(NUP,NI0)  00001179 


COMMON  /STORE/  AE , A A , BB , BL , CL , T V , T8 V 
COMMON  /8T0R2/  IMX 

common  /Insub/  TL («) ,MT,TOW,ON,STARTR,OELTR»RN(a) ,8C0N,CElhTH(3) 
COMMON/AZM/  AZMUTH(460)  ,NA,'ELEVAT,PRF,KEEP 
COMMON  /AOATA/  10 A V , IMOUR r IMIN , ISEC , NTP, NSF , NOD , NRC 
C0MM0N/w0RK/W(5lfl);v(5U),V8(51«),SV(5ia),V8(5l4),VJ(5ia),UI(5ia) 
IVSTtBia) ,Hl {5l«),H2(5ia) ,NCL 
common  /AZ2/8AZ,CAZ,nAZ, I8CANF,nEL 


•00001160 


IEM  is  no, of  events  in  c radial, 
initialize  and  generate  hc  array 


1EM«XCVNT(1) 

LM«5*3oITr 

LMMiLMal 

IDX»LM^1 

NCLMiNCL«l 

L0MX«(NUMAX.2)/LM 

IF  (LDB.GT.LOMX)  L0B»LDMX 

numP«2»LM*L08 


zero  arrays 


DO  11  I«1,KMAX 
DO  11  J*1,JMAX 
IPL0(J,I)«0 
DO  21  IB1,IEMAX 
IPTC(n»0 
DO  21  KM1,KMAX 
TC(K,I)»0 
IPCNT(K,I)«0 
DO  21  J«1,JMAX 
IPCl  (J,KiI)»0 
IPC2(J,K.I)»0 
1PC3(J,K,I)»0 


00001 IBl 
00001162 
00001163 
00001164 
00001165 
,00001166 
00001167 
00001166 
•00001169 
00001190 
00001191 
00001192 
00001193 
00001194 
00001195 
00001  196 
00001197 
00001196 
00001199 
00001200 
00001201 
00001202 
00001203 
00001204 
00001205 
00001206 
00001207 
00001206 
00001209 
00001210 
00001211 
00001212 
00001213 
00001214 
00001215 
00001216 
00001217 
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T17 

8B1 

921 


912 


714 


933 


X 2X,4H(KM),SX,lH{0BZ)*6X,llHeDBZ,KMi»*2),l)(i7M(M/8EC),lX« 

x7H(H/5EC),1X,11H(M/9EC)**2)) 

TZ«TZ/I0e05 

WRITEC6»7l7)X,AvZ,TZ»ve, VN;vER.OEL 

FQ«MAT(If»kX»F5,l»7X,F9,l,fX,>fc.l,2X,Fl,,l,2X,Fl0.l,lX,ll5,S) 
W4ITE(6,«S1) 

F09MATtlX,«(*— ) 

CONTINUE 
ID1««999 

MRITE(2)XBAR,VBA9,ID1 
CALL  FA6E 
W91TE(»,932) 

FORMATUH  ,*  PEAK  OETECTEO  CELL  ATTRIBUTES*) 

WRITEC6,714) 

FORM AT (41 Xf7HAVERAQE*3X«7HA VCR age/sox, BHLOC AT  ION, SX|6HR ADI AL,2X, 
X10HTAN6EntIAL#1X,10HTAN6CNTIAL/BX,12HPEFLECTIVITY,3X,4HAREA,Sx, 
X4HEA8T,2X,9HN0RTH,3X,SH8HEAR,4X,SH8HEAR,6X,SHSHEAR/2X,2HI0,5X, 
X5H(D8Z) ,SX,TH(KM**2),2X,4H(KM) ,SX,  4H (KM) , BH (M/s/KM) , 3X , 
X§H(M/8/KM) ,3X,BM(M/S/KM) J 
DO  933  Nil.NCEU 

IF(UPC1#N),LE.O,,OR.UP(2,N),EO.O.)60  TO  933 
UP(I,N)«UP(3,N)/UP(2,n1/10E08 
UP(4,N)BIJP(4,N)/UP(2,N)/IQC02 
UP(2,N)RUP(2,N)/UP(1,N) 

UP(5,N)PUP(3,N)/UP(1,N)*10E02 
UP(6,N)»UP(fc,N)/UP(l,N) 

UP(1,N)RUP(1,N)*DELR/1.E06 

wrITC(6,718)N,UP(2,N) ,UP(l,N),UP(3,N),uP(4,Ni,UP(5,N) ,UP(fc,N) , 
XUP(7,N) 

718  FORMAT(l4,SX,F5,l,feX,F6,l,lx;Ffe.l,lX;F6.1,lX,F7.2,4X,Fb,2,4X,F7,2) 
CONTINUE 
CALL  PAGE 
wRITE(b,942) 

FORMATdH  ,*  TANGENTIAL  SHEAR  MAXIMA  ATTRIBUTES*) 


942 

715 


719 

943 

950 


wrITE(6,715) 

F0RMAT(25X,8HL0CATI0N/TX,SHSHEAR,4X,4HAREA,3X,4HCA8T,2X,5HN0RTH/ 

X2HI0, 

XlX,8H(M/9/KM) ,1X,7H(KM**2) ,2X,4H(KM) ,3X,4H(KM)) 

DO  943  N«1,NVCEL 

IF(UV(1,N) ,LE,0,,0R,UV(2,N) ,EO,0,)GO  TO  943 
UV(4,N)BUV(4,N)/UV(2,N) 

UV(S,N)«UV(3#N)/IIV(2,N) 

UV(2,N)«UV(2,N)/UV(1,N) 

l)V(l,N)PllV(l,N)*OEl,R/l.EOb 

WRITE(b,7l9)N,UV(f  ,N)  ,UV(1,N)  ,IIV(3,N),IJV(4,N) 
F0RMAT(l4,lX,F7.1,3X.Fb,l»lX,Fb,lilX*Fb,l) 

continue 

WRITECk, 950)100 
FORMATdHO,  IOHTOTAL  IDD»,I6) 

I8CANF*0 

RETURN 

End 


00001 

00001 

00001 

00001 

00001 

00001 

00001 

00001 

00001 

00001 

00001 

00001 

00001 

00001 

00001 

00001 

00001 

00001 

00001 

00001 

00001 

00001 

00001 

00001 

00001 

00001 

00001 

00001 

00001 

00001 

00001 

00001 

00001 

00001 

00001 

00001 

00001 

00001 

00001 

00001 

00001 

00001 

00001 

00001 

00001 

00001 

00001 

0000} 

00001 

00001 

00001 

00001 


111 

112 

113 

114 

115 

116 
117 
116 
119 
ICO 
ICl 
122 

123 

124 

125 

126 

127 

128 

129 

130 

131 

132 

133 

134 

135 

136 

137 
13B 

139 

140 

141 

142 

143 

144 

145 

146 
14? 
146 

149 

150 

151 

152 

153 

154 

155 

156 

157 
15B 

159 

160 
161 
162 
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671  DEL*)a8C0N*CELWTH(NTP^l) 

WI»ITE(2)X,V,IPU 
CALL  PAGE 
wPITEC6,87?) 

872  POPMATflH  FIXED  CONTOUR  ATTRIBUTES*) 

WRITE(8,712) 

712  FORHAT(28X*7HAVERA6E«SX,8HLOCATION,ax,9HTOTAL»3Xf 7HAVERA6E/ 
X^X,9HTHRE8N0LD,SV,aHAREA,2Xi 12HREFLCCTIVITY»1X,4HEA8T»2X, 
X8HN0RTH,2X»8HPREClP,ax,6HPRCCIP/2X«2HI0»aXf 

X5H(DBZ) ,ttX»7H(KMft*2) ,4X,SH(DBZ) ,SX,4H(KM) ,1X,aH(KH) , 1 X • 9H (T0N8/HR) 
X1X,7H(MH/HR) ) 

DO  931  KPl.NFC 
JOOiFDO(K) 

00  931  JtliJDO 
IEOGEpO 

IF  (ATR(IAT,J|K).E0,0,0)  GO  TO  931 
IF  (ATRdAT,  J,K)  ,LT,0,)  IE0GE>1 
lOlPJ 

A8ARpDCLR*ATR(1, J,K) 

IF(ABAR.LE.0.)60  to  931 
ZBARpATR(2i J»K)*0ELR/ABAR 
X8ARpATR(3i J»K)*0ELR/ABAR/ZBAR 

VBARpATRC4i j»k)*oelr/abar/zbar 
ABAR9ARAR/1000**2 

IF  tK.GTtl) write (fc, 720) 101 iTL fK) ,ABAR,ZBAR,X8AR,YBAR 
720  F0RHAT(1X,IS,5X,I2.'4X,F9,2,4X,F5,1.«X,2F6,1) 

IF  (K.GT.n  GO  TO  931 

TP»ECp08I(J)*0ELR 

AVPRECaTPRCC/(ABAR*1000) 

TPRECpTPREC/IOEOS 

wrItE(6»716)  I01,TLfK)  , ABAR.'ZBAR.XBAR,  YBAR,TRREC,AVPREC 
716  format (1X,I3»5X,  12. ''9.2, «X,FB,1»4F,F6.1»1X,F6.1.1X»F7,2,2X, 
XF6,2/) 

YBAR«YBAR/(1,84*10E03)*4,0 
X8ARpXBAR/(3.8A*10E03)*4.0 
WRITE(f)XBAR,YBAR,IDl 
931  CONTINUE 

00  921  IPUNNE 

IFfZM(NZP,I,J),CO,0)GO  TO  921 
TZpDELR*ZH(11,I, J) 

AVZrZH(11,I, J)/ZH(12,I, J) 

C COMPUTE  AVO  WINO  8PEE0  AND  OIR, 

DEL»ZHf5,I, J)*ZH(6,I, J).ZHf7,l, J)*ZH(7;i,  J) 

IF  (OEL.EQ.O.)  go  to  921 

VNP(ZM(4,1, J)*ZH{5,I, J).ZMC3,I, J)*ZH(7^I, J))/DEL 
VEp(ZH(6,I, J)*ZH(3,I, J)-ZM{7,I, J)*IM(4,I, J))/0EL 

VERP(ZH(2,I, J)*VN*VN*ZH(6.I» J)*VE*VE*ZM(5,I, J)*2.0*VN*VE*ZMC7,I, J) 
1«2.0*VN*ZH(4,I, J) *2.0* VE*ZH (3,1, J))/ZH (10,1, J) 

VF»ZH(5,I, J)/ZH( 10,1, J)-(ZH (8,1, J)/ZH( 10,1, J)) 

VC»ZM(6,  I, J)/ZH(10, I, J)-(ZH(9,1, J)/ZH(10,1, J))**2 
WRITE (6,713) 

713  FORMATdHO,  10X,7HAVERA6E,8X,5HT0TAL.9X,'7HAVERA6E,1X,7HAVCRA0E,2X.' 
X8HvEL0CITY/1X,6MHCight,1X,12HREFLECTIVITY,  1X,12HREFLECTIVITY,4X, 
X1H(I,7X,  lHV,BX,8MVARIANCE,iaX,3H0EL/ 


000010S6 

000010S7 

00001058 

00001059 

00001060 

00001061 

00001062 

00001063 

00001064 

00001065 

00001066 

00001067 

00001068 

00001069 

00001070 

00001071 

00001072 

00001073 

00001074 

00001075 

00001076 

00001077 

00001078 

00001079 

00001080 

00001061 

00001082 

00001083 

00001084 

00001065 

00001086 

00001087 

00001068 

00001089 

00001090 

00001091 

00001092 

00001093 

00001094 

00001095 

00001096 

00001097 

00001098 

00001099 

00001100 

00001101 

00001102 

00001103 

00001104 

00001105 

00001106 

00001107 

00001108 

00001109 

00001110 
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NAMEal2 

00001328 

WRITe(6,9909)NAME,MPCEL 

00001329 

c 

00001330 

c 

ASSOCIATE  CELLS  ON  B RaDXAL*  TOP  DOWN 

00001331 

c 

00001332 

193 

mpkro 

00001313 

NAMEtll 

00001334 ‘ 

write (6, 9909) NAME, NPK 

00001335 

9909 

F0RMAT(5)(,I2,10K,IS) 

00001336 

IPCNA.EQ.nSO  TO  361 

00001317 

TATM»0, 

0000133a 

iHBrlPCl (IPE,KC,IE)61 

00001339 

IHDfIPC2fIPE,KC,IE) 

00001340 

DO  194  I»IHB,IHD 

00001341 

IFfHSd)  , BO, ••999)60  TO  194 

00001342 

If  (IAB8(MB(I)),6T,TC(KC,IE)4L0B)60  TO  934 

00001343' 

194 

CONTINUE 

00001344 

If  (JE2«E0«0)  60  TO  361 

00001349 

DO  261  je«JEl,JE2 

00001346 

If  CIBt2, JE, 1) .LT.IPCl fIPE,KC,lE))  SO  TO  261 

00001347 

If  (IB(1,JE,1),6T,IPC2(IPE,KC,IE))  GO  TO  361 

00001348 

C 

00001349 

C 

JE  EVENT  ON  B RADIAL  IS  ASSOCIATED 

00001390 

C 

00001391 

271 

IPBbIPTBCJE) 

00001391 

00  291  L0Rl,lPB 

00001393 

t(SaIPR«LB«l 

00001394 

NPl*lPPNT(t<8,  JE) 

00001399 

DO  281  JPE«l,NPl 

00001396 

If  (IPR2(JPE,K8, JE) .LT.IPCl (IPE,KC,IE))  GO  TO  281 

00001397 

If  (IP8l(JPE,KB,JE),6T,IPC2(IPE,KC,IE))  SO  TO  361 

00001398 

LPCEL»IPB3CjPE,K8, JE) 

00001399 

WRITE(6,2729)  JPE,IPE,KB,KC,  JE,IE,LPCEL.’HPf 

00001360 

2729 

f09NAT(2X,8I10) 

00001361 

If  (LPCEL,Ea,0)SO  TO  281 

00001362 

TATMiAMAXl (TATM,TATR(1,LPCEL) ) 

00001363 

If  (TATM,E0,TATRC1,LPCEU))NPK«LPCEL 

00001364 

291 

continue 

00001369 

291 

CONTINUE 

00001366 

261 

CONTINUE 

00001367 

If (MPK,Ea,0)60  TO  361 

00001368 

If  (AB8(TATR(1,MPK))  ,6T,TC(»<C,IE)*L0B)MPK*.MPx 

00001369 

SO  TO  361 

00001370 

934 

MPKA-CNID+l) 

00001371 

C 

00001372 

C 

have  B compare  WITHIN  RAN6E 

00001373 

c 

00001374 

361 

continue 

00001379 

NAHEa2l 

00001376 

WRITE(6,9909)NAME,MPK 

00001377 

If  (MPK,E0,0,An0,NPK,EQ.0)S0  to  631 

00001378 

C 

00001379 

C 

MPKwO,ANO,NPK«0  «»  NO  COMPARE 

00001380 

c 

HPRbO,ANO^NPK,NE<0  * no  8 COMPARE 

00001381 

00001382^ 

c 

NPKiO.ANO.MPK.NC,n  • B COMPARE 

!F  tl4B8CUtI»n).l.f,TC(K,XE)l  Oo  TO  12J 

00001273 

SO  TO  I3l 

0000127G: 

c 

0000127S 

c 

START  RANGE  FOR  SEGMENT  (CONTOUR) 

00001276 

c 

00001277 

121 

IPCNT(K,XE)«XPCNT(t<,IE)*l 

00001276 

IPE»IPCNTCK,IE) 

00001279' 

IPCI(XPE,K,IE)«I-1 

000012SO 

IPLO(XPE,K)«XP»(.l 

000012S1 

131 

continue 

00001282 

60  TO  161 

00001283 

C 

00001284 

c 

END  RANGE  FOR  SEGMENT 

0000128S 

c 

00001286 

141 

no  151  KL«K,IPT 

00001287 . 

IF  (U(I-1),EQ.»999)  SO  TO  161 

000012881 

IF  (IAB8(U(X-n)  ,LE,TC(KL»IE))  GO  TO  161 

00001289 

IPE«IPCNT(KL»IE) 

00001290 

XPC2(IPE,KL,  lE)m«l 

00001291 

151 

continue 

00001292 

161 

CONTINUE 

00001293 

00  181  !<«1,IPT 

00001294, 

XPE»IPCNT(K,IE) 

00001299 

00  181  miflPE 

00001296 

WRlTE(6#17niE,I,K;iPCl(l,K,IE),IPC2(I,‘K,lE) 

,IPCNT(K,IE),IPL0(I, 

K)00001297 

liTC(K,XE) 

00001298 

ITl 

FORMATdH  ,3X3,5110) 

00001299 

181 

CONTINUE 

00001300 

C 

00001301 

C 

associate  CELLS  LOOP  ON  THRESHOLD 

HIGHEST  TO  LOWEST 

00001302 : 

C 

00001303 

940 

00  941  LC«l,IPT 

00001304 

KCaIPTIpLC#! 

00001309 

NPC«IPCNT(KC»IE) 

00001306 

C 

LOOP  ON  SEGMENTS 

00001307 

00  941  IPE»1»NPC 

00001308 

KbKC^I 

00001309 

NPKaO 

00001310 

TATM«0, 

00001311 

IF  (K.GT.IPT)  GO  TO  193 

00001312 

LPE«1PCNT(M,IE) 

00001313 

192 

00  191  Lil»LPE 

00001314 

IF  (IPC2(L,K,IE) .LT.IPCl (IPE,KC, IE))  GO  TO  191 

00001315 

IFdPCl  (L,K,IE)  ,GT,IPC2(IPE,KC.1E))GP  TO  193 

00001316 

NPCELmIPC3(L»K,IE) 

00001317 

IF  (NPCEL.EQ.O)  GO  TO  932 

00001318 

TATM«AMAyi(TATM,TATR(l,’NPCEL)) 

00001319 

IF(T4TM,E0,TATR(1,NPCEL))NPK«NPcEL 

00001320 

C 

00001321 

C 

NPCEL  IS  FOR  NEXT  HIGHER  (ENCLOSED) 

THRESHOLD  ON  C RADI ALOOOOl 322 

C 

00001323 

231 

IF  (AB8(TATR(1, NPCEL)). CT,(TC(KC.IE)*LOB  )) 

60  TO  932 

00001324 

191 

continue 

00001329 

60  TO  193 

00001326 

932 

NPKa«l 

00001327 

89 


c 

c 


c 

c 

c 

301 


392 

359 

391 


aoi 


an 

419 


HIGHEST  THIS  940X41 


Xf tHPK,CQ,0,4N0,NPK,LT,05G0  TO  931 
IP(MPK,NE,0)00  TO  421 


NO  PRIOR  R4DI4L  EOR  C0MP4RXS0H,  XNCREHEHT  NPClL 


359 


C 

C 

c 

366 


341 

391 


NPCEU«HPK 
IP(N4,EQ,nG0  TO 
00  352  X«IHB,IHD 
XPfH8(l),EQ, "999)60  TO  392 
If  (X4B8(HBCI)) ,GE,TC(KC,XE))GO  TO  931 

continue 

IN0X»T4TR(1 ,NPCEL)-TC(KC.IE)"1 

If  {INOX.GE.LOS)  GO  TO  931 

IPC3(IPE,KC,IE)«NPCEL 

INP1*IN0X*LH 

INX»10X  + IN0Xi»LM 

If(N4,EQ.nG0  TO  419 

I8T*IPC1(IPE,KC,IE)*1 

18P«XPC2(IPE,KC,IEi 

DO  411  IPI8T,I8P 

R»8C0N*(Fl04T(I"1)",9)*CELWTH(NtP*1) 
T4TR{Xnt1,NPCEL)«T4TR(INt1,‘nPCEl)t04Z*R 
T4TR(1nT2,NPCEL)pT4TR(INt2,NPCEl)+DAZ*RpU(I) 
T4TRf Int3,npCEL)«T4TR(TN*1iNPCEl)*0AZ*84Z4R*I 


If 

If 


(ITV.ne.I)  go  to  411 

(vCI),E(5,"999,0R,V(Ip1),E(J,"999)  60  TO  40J 


If  CVSd)  ,E0."999)  GO  TO  411 
T4TRCInT6,NPCEL)«T4TR(IN46,NPCEl)*R*VS(I) 


CONTINUE 


0000136S 

OOOOlSli 

00001309 

00001304 

00001387, 

00001380 

00001389 

00001390 

00001391] 

00001390 

00001393 

00001398 

00001399 

00001394' 

0000139TJ 

mm  mi 


00001398 

00001399 

00001400 

00001401' 

00001400 

00001409 

00001400; 

00001409 

000014041 

00001407 


10) 

N4HEB31 

wRITE(6,9909)N4ME, INOX 


U(I) 

00001408 

U(I) 

00001409 

00001410 

00001411 

V(I"1)) 

00001410 

00001413 

00001414 

IAB8(VS(X)))) 

00001419 

0000)414 

) 

00001417 

tTATR(lNX,NPCEL);-l 

1.00001418 

00001419 

000014201 

DO  365  L«1,LPE 


LPCEL»IPC3(L»f»XE) 
If (LPCEI.E0,0)G0  TO 


931 


If  fNPCEL.EO.LPCEDOn  TO 


(KZ, 

NPCEl) 

»KZ 

■l.NUMP) 

ITH 

NPCEL 

AT 

THI8 

level 

IfE, 

9Ci 

IE) 

)60 

TO 

365 

IRE, 

KC. 

IE) 

)Gn 

TO 

931 

)on 

TO 

365 

369 

,1E) 

"1 

000014211 

00001400 

00001421 

00001420 


00001409' 


365 


INX«IDXTINDX*LM 
If fINX,6T,NUMP)G0  TO 
T4TR(lNX»LfCEL)»0, 
T4TB(2*INDX*LH,LfCEL)*NPCEL 


00001424 

00001427 

00001420 

00001429 

00001430] 

00001431 

00001430 

00001433 

00001434' 

00001439 

00001434 

00001437 
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o o o 


XNOHLOB«XNOX 
00  «83  XlliXNO 
DO  «BS  Jll*LH 
XNaUJ*(LOB«X)*LH 
XMaiaJa(XND»X)*LH 
4B3  TATR(XN,LPCEL)aTATR(IM;L?CELj 
TNOtXNDX*LH*l 
4B2  00  4B4  XlXNS»XND 

4B4  TATR(X,LPCEL)«0. 

4BB  XNaO 

XPC3(XPEiKCiXE)«LPCEL 

NAMEaBl 

write (6*4904) NAME fXNDX 
60  TO  Si2 

4B5  DO  4B6  iRiiNlO 

XP(IACT(I) .10.0)60  TO  4B7 
466  CONTINUE 

WRITE(6*644) 

GO  TO  431 
487  LPCElal 
XACTCX)*! 

TATB(1*LPCEL)wTC(KC,IE)*1 

TATR(NUMP,LPCEL)aXC(NPA,ie,l) 

60  TO  468 

422  LPCEI«XABS(MPK) 

1F(LPCE1.6T.NX0)B0  TO  431 
DO  441  JEaJEl*JE2 

IF  CIB(2,JE,1).LT.IPCI(IPE*KC»IE))  GO  TO  441 
IF  (IB(1,JE,1).GT,I6C2(IPE,KC,IE))  GO  TO  632 
IFB4IPTB(JE) 

DO  471  LB«1,1P8 
KBaIPB»LB«l 
MFB«IPBNT(KB* JE) 

DO  461  JPErIiMPB 

IF  (IP82(JPE»KB,JE),LT,IPCI(IPE,KC,IE))  GO  TO  461 
IF  (IPBI(JPE»KB,JE).GT.IPC2(IPE,KC,IE))  go  to  632 
NPCELwlPB3(JPE,K8, JE) 

IF  (NPCEL.LE.O)  go  to  461 
IFCLPCEL.EQ.NPCEDGO  TO  461 
XF(TB(KB,JE).GT.TC(KC,IE))80  to  461 

COMBINE  AT  T8»TC  LEVEL  AND  BELOW 

902  IndX«TATR(1,NPCEL)»T8(K8, JE) 

NAHEbTI 

WRITE(6,9404)NAME.INDX 
IF(InDX.6E.L0B)60  to  461 
XMOXaTATR(l.LPCEL)«TB(KB, JE) 

NAHEaBl 

WRITE(6|4909)NAMEf IMOX 
IF(ImDX.LT.LDB)GO  to  661 

891  DO  692  Jai.NUMR 

892  TATR(J,NPCEL)»0. 

IACT(NFCEL)aO 

NAMEalOl 


00001443 
00001444  j 
00001449 
00001446  : 
00001447  ; 
00001498 
00001444 
00001900 
00001901  . 
00001902 
00001903 
00001904  : 
00001909 
00001906 
00001907  . 
00001906 
00001904 
00001910 
00001911 
00001912 
00001913 
00001914 
00001919 
00001516  ; 
00001517 
00001518 
00001514 
00001920  ; 
00001921  ^ 
00001522 
00001923 
00001524  ; 
00001525 
00001526  : 
00001527  j 
00001526 
00001524  ' 
00001530 
00001531 
00001532  i 
00001533 
00001934  : 
00001939  ' 
00001936 
00001537 
00001938 
00001534 
00001940 
00001941  ^ 
00001542 
00001S4S 
00001544 
00001549 
00001946 
00001947 
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write  (b»  ^90^5  NAHf,igpcEL 
60  TO  a61 
861  IN0«ZNDX«LI)B 
DO  691  N«1,IND 
LDll^(Lt)B*M)ft|.H 
WORl«(ZNDX»M)ftLM 
DO  691  Is1,LM 
IF  (IteC.LH)  60  TO  661 
IF  (I,6E,8)  80  TO  871 

TATRCL06I,LPCEL)"T4TRCN0*I,NFCEl)*T*TR(L0^I,LFCEL) 

TATR(NO^Z,NPCEL)bO,' 

60  TO  691 

871  TATR(LD6I,LFCEL)»AMAXI  f TATR (nO*I , NPCELJ  , TATR fl.D^I,LPCEL) ) 
TATR(NO^I,NPCEL)RO.‘ 

80  TO  891 

861  TATR(NO,I)iLPCEL 
891  continue 
461  continue 
471  continue 
441  continue 
632  IFtNPK.LE.Ojeo  TO  931 

npcelrlpcel 

60  TO  366 
C 

C UNASSQCZATEO 

C 

631  IF(NA.EQ.1)60  TO  639 
00  641  I«IHB,IHD 
IF  (HBCn  ,EQ,«999)  60  TO  641 
IF  {IAB3(HB{I))  ,GE.TC(KC,IE))60  TO  931 

641  CONTINUE 

639  00  642  Jai,NlD 

IF(1ACT(J),EQ.0)6D  to  643 

642  continue 
WRITE(6f644) 

644  F0RMAT(5X,*  TOO  MANY  CELLS*) 

GO  TO  931 

643  NPCEL«J 
1ACT(J)P1 

661  IPC3(lPE,KC,ZE)aNPCEL 
IPKilPLOtlPE.KC) 

IRfIPCRNG(IPK,IE) 

IN1«LN*1 

IN»(LDB«n*LM*lNl 
00  671  1«IN1,IN 
TATR(1,NPCEL)*0,0 
671  continue 

591  TATRri»NPCEL)*IAB8fUflR)) 

TATR(NUMP,nPCEL)»IC(NPA,IE, 1) 

IFINA.EG.DGO  to  939 

I8T»IPC1(IPE,KC,IE)*1 

ISP*IPC2(1PE»KC,1E) 

00  621  l«18Tfl8P 

RrSCON* (FLOAT! 1-1) •,5)*CELWTH(NTP*1) 
TATR(2,NPCEL)«0AZ*«*TATR(2,NPCEl) 


0000 

0000 

0000 

0000 

0000 

0000 

0000 

0000 

0000 

0000 

0000 

0000 

0000 

0000 

0000 

0000 

0000 

0000 

0000 

0000 

0000 

0000 

0000 

0000 

0000 

0000 

0000 

0000 

0000 

0000 

0000 

0000 

0000 

0000 

0000 

0000 

0000 

0000 

0000 

0000 

0000 

0000 

0000 

0000 

0000 

0000 

0000 

0000 

0000 

0000 

0000 

0000 

0000 

0000 

0000 


92 


r 


TA79t3,NPCEL)»0A2*R*U(I)4T*TR(3,NPCEL) 

0000160] 

TAT9(4,NPCEL)»0AZ*6Af*»*R*l'(I)4TAT9(4,wPCEL) 

0000160A 

TATR(5,NPCEL)»DAZ*CA7*»*9*U(n*TATR(5,NPCEL) 

00001605 

IF  (ITr.Ng.n  GO  TO  621 

00001606 

IF  CV(1),EQ,»999,0»,VC1«1),EO,.999)  GO  TO  601 

00001607 

TATR(6,NPCEU)«0AZ*R*(V(n-V(I-n)*TAT9f6,NPCEL) 

00001608 

601 

IF  (V8(l) ,E0,.999)  60  TO  621 

00001609 

TATR(7,NPCEL)«DAZ*9*V8tn*TATR(T»NFCELi 

00001610 

TATR(8,NPCEL)«AH4xi (tATR(8,NPCEl)*FL0ATCIA83(V8(I)))) 

00001611 

621 

CONTINUE 

00001612 

6S9 

TATPCI0X,NPCEL)«NA 

0000161S 

IFfI8T,E0,2,0R,ISP,EQ,lHX)TATR(l0X,NPCrL)»»TATR(IDX,NPCEI.) 

00001614 

WRITe(6,l07nNPCeLt  (TATR(KZ,nPCEL5  |KZ«1,NU«P1 

00001615 

931 

COMTINUE 

00001616 

941 

continue 

00001617 

951 

continue 

00001618 

IF(NA,EQ.nG0  TO  1051 

00001619 

C 

00001620 

C 

end  of  association  loops 

00001621 

C 

00001622 

ID2p1+(L08»1)*LM 

00001623 

Lnx»l^LO0*LM 

00001624 

00  991  iPliNlO 

00001628 

IFflACT(I) ,EQ,0)G0  TO  991 

00001626 

IFfTAT»(l,n,GT,0,,AN0,TATR(2,I),GT,0,iG0  TO  961 

00001627 

GO  TO  991 

00001628 

961 

IF(TAT«(LDX,I).LE,O,,0R.TATRfin2*l,I),LE,0,)GO  TO  991 

00001629 

IF  (TATR(LOX,n,LT,NA,nR,ISCANF,EQ,nGn  TO  971 

00001630 

GO  TO  991 

00001631 

971 

DO  961  J«1,LNN 

00001632 

LlPf  J,NCELL)»TATR(I02*J,n 

00001633 

961 

CONTINUE 

00001634 

UPfLH»NCELL)PTAT9(NUMP,n 

00001635 

NAMEPlOl 

00001636 

WRITE(6,9909)NAME, I 

00001637 

write  (6»  1071)  I,  (TATRfK,  I)  ,t<»l  ,NUHP) 

00001638 

wRITE(6,9910)NCELL,CUPfK,NCELL) .KbI ,LH) 

00001639 

9910 

F0RMATtlX,12,12X,8Fl5.2) 

00001640 

NCELLPNCCLL+l 

00001641 

DO  982  JB1,NUNP 

00001642 

962 

TATRIJ,  n«0. 

00001643 

1ACT(1)pO 

00001644 

991 

CONTINUE 

00001645 

1031 

00  1041  IbI,1EMAX 

00001646 

IPTB(I)bIPTC(I) 

00001647 

00  1041  Kb1,KMAX 

00001648 

T8(K,I)BTC(K,n 

00001649 

IP8NT(K,I)BirCi>IT(K;i) 

00001650 

00  1041  jBljJMAX 

00001651 

IPBl  CJ,K,I)bIPCI  (J,'',1) 

00001652 

IP82(J.K>I)BlPC2(J,K,n 

00001653 

IPR5(J,K,nBlPC3(J,K,l) 

00001654 

1041 

CONTINUE 

00001655 

DO  1 Ib2,NCLM 

00001656 

MMa«999 

00001657 

I 

I 

i 


If  (U(I-l)  ,NE,.9995  MH«IAB8(U(Ul)) 

If  (U(I) ,NE,-999)  MHbMAXOCMH, lA88(U(n ) ) 

IF  tf''E,«999)  MH«MAX0(MH,1ABS(U(I«»1))) 

I H0(n«MH 

Nal 

write  (6,1091)  U 
1061  FORMAT  (16) 

DO  1081  1*1, NID 

IF(IACTCI) ,CQ,0)GO  TO  1081 

WRITE(6, lOTl)  I, (T6TR(K,I),K«1,NUMP) 

1071  FORMATdX,  ia,SX,9Fl3,2,/(19X,8Fl3,2)) 

1081  continue 

Nag 

WRITE(6, 1061)N 
IC82  no  llOl  IE»1,IEM 
IPT»IPTB(IE) 

DO  1101  Kal,lPT 
IPE«IPRNT(K,IE) 
on  1101  Ial,IPE 
1T*TR«0 
TATRX*0, 

IPX»IP83(I,K, I£) 

IF(IP*,GT,0)TATRXaTATR(l ,IPX) 

IF(IPX,GT,0)ITATR«TATR(IOX,IPX) 

wRITE(6,1O91)I,K,IE,IP0i (I.k.IE) ,IP82(T,K,IE),IPB3(I,K,IE),TB(K, 
HE)  ,TAT»X,|T6T» 

1091  FORMATdH  , 313,  al8 , E I S , 3 , 1 8 ) 

1101  continue 

RETURN 

Eno 


0000165 

0000165 

0000166 

0000166 

0000166 

0000166 

0000166 

0000166 

0000166 

0000166 

0000166 

0000166 

000016? 

0000167 

0000167 

0000167 

0000167 

0000167 

0000167 

0000167 

0000167 

0000167 

0000168 

0000168 

0000168 

0000168 

0000168 

0000168 

0000166 

0000168 
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I 
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MICROCOPY  RESOLUTION  TEST  CMA5?T 

NATIONAL  BURtAlJ  Of  STANOAROS 


o n o oo 


« 


SUBROUTINE  PRN}  OOOOUBi 

****ft***************;*;***;****«******;*****;**********ft********;*oooous| 

PRINTS  OUT  UNPACKED  DATA  (INTCGCR  P0RNAT)»  00001690 

version  1.0  LEVEL  7609E0  OOOOUff 

JHW  CDC  6600  APSL  PE09S  OOOOUfi 

*************tt******;*******;****«*«**;******«**«**ft******ft***;*;*ooooi69| 

C0NM0N/A1026/HVP(S,'l0E4)  00001694 

C0NM0N/AZM/AZHUTH(460)«NA.ELEVAT«PRF.KEEP  0000169S 

CONMON/AOATA/IDAyilHOUR^ZMIN.ZSEC.NTP.NSF.NODfNRC  00001694i 

C ..•.•.•.•..•.•0000l69f 

CALL  PAGE  00001694 

write  (6»  100)  ID  at.  IHOUR*  ININ,  ISeCfNTP.'NSFf  NOD.  NRCfPRFfAZMUTH(NA)f  0000169i 
X ELEVAT  00001704 

100  FORNATCSSH  unpacked  radar  data  (INTEGER  VALUE8)/2Xf  00001701 

XAOAT  hr  MN  sc  CELLWIOTH  subframe  dump  FREO  «CELLS  PRF  AZIMUTH0000170I 
X ELEVATI0N*/I9.Sn.9x;i2»7x;i2.BX,l2,I7.F6.0»Fl0.1.SX,F9.1//4X>  00001701 

X1H1|2X,B(14H  mean  VAR  PWR)/)  00001704 

KKVKEEP  0000170S 

NNMKEEP67  00001704 

DO  10  Nil, 32  00001707< 

WRITE(6,101)KK,  ((MVPd,  J),Iil,3) , JaKK,NN)  00001704; 

101  FQRMAT(I9,2XiS(I6,21«))  00001709 

KKiKK^S  000017109 

NNiNN^S  000017111 

10  continue  ooooiTii: 

RETURN  00001711 

END  00001714 


o o o no 


10 


100 


SUBROUTXNK  R»NS(M00E,W) 

PRINTS  B8CAN  HAPS  OP  COMPUTED  AND  CODED  DBZ  AND  VELOCITY. 
VERSION  1,0  LEVEL  7B11I9 

JHW  CDC  B600  APBL 

INTEBER  TL»8TARTR»w 

logical  PRlNTl,PRlNTE»PRINTJ,PRlNT«fCONTRZ,CONTRV 

00001715 

)*t000017lB 

00001717 

0000171S 

00001719 

»*«0000t7t0 

000017!! 

000017!! 

DIMENSION  IC(64),Wfl) 

000017E4 

COMMON/PARM/PRINT}, PRINT!. PRINTS, PRINTS. ICODES (SB). At. SI.AE. BE, 

00001715 

contre,contrv.npile.nrec.numr 

000017EB 

C0MM0N/A10Z4/  MVP(S.10!4) 

000017E7 

COmhON/AZH/  AZMUTHf4B0),NA.ELEVAT.PRP.KEEP 

000017E6 

commqn/adata/idat.ihour.ihin.isec.ntp.nsp.ndd. 

NRC 

000017E9 

.IsAnnAt  V Vn 

NC0L"(NRC41)*!9B 

00001791 

f 

f 

INT4NC0L/B4 

0000179! 

IP(M00E.E0,1)CALL  cohpz 

00001799 

!■! 

00001794 

L41 

00001795 

i 

iViO 

00001796 

IP«0 

00001797 

DO  10  N«t,XNT 

00001798 

1 

L»L^1 

00001799 

( 

IV4XV»N(L) 

00001740 

1 

[ 

IP4lP*l 

0000174} 

/ 

IVBIV/IP 

0000174! 

1 

ir«Al*IV  ♦ B1 

00001749 

h 

XP(IT.QT.S6)IY«SB 

00001744 

IP(IT,LE.0)IV»1 

00001745 

! 

xccn^icoDEsdv) 

00001746 

1 

I4141 

00001747 

IP(L.LT,NC0L)00  to  3 

0000174S 

i 

WRITECfc,10O)AZMUTH(NA),ELEVAT.X0AY,lH0uR,IMXN, 

ISEC.XC.PRP 

00001749 

F0RMAT(1X,PS.1,PB.1.I4,1X.2Z!.XS.9X.64A1.SX.P7 

• 1) 

00001750 

ii 

RETURN 

00001751 

i- 

END 

0000175! 

f; 
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i 


o n n n 


SUBROUTINE  PA6E 


PRINTS  PAGE  HEADER  AND  KEEPS  TRACK  OP  LINE  COUNT  00001755 

VERSION  t.O  level  711122  0000175S 

***********«******;*******************;**ft**;********ft***«****;*;*ooooi7S7 
integer  ICOOEiIRUN.NPAGE  0C00175S 

REAL  TITLE(6)  0000175R 

[%  COMMON  /HEAD/  TITLE, ICODEf VERS. LEVELfDATi  IRUN, nPAGE, NLOG  00001760 

COMMON/LINUM/LINE  00001761 

LINEtA  00001763 

[ NPAGEiNPAGEtl  00001764 

write  (6i2030)  ICOOE.IRUN, title, VCRS, level, DAT,  NPAGE  00001765 

2030  P0RMAT(wl*,lS,I6,5K,6A6,*  VERSION  *,PS.l,*  (*, 16,*)*, IIK,  00001766 

, X A10,10X,*PA6C  *,'I3/1X,'127(***))  00001767 

RETURN  00001766 

END  0000176R 


; 
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subroutine  LXNES(N),RETURNS(A)  00001770 

C *********************«**********************************t*****ftft;*00001TTl 
C VERSION  1.0  lEVCt  7S09I1  00001772 

c **********«*************************«********«*«**«**t********;*;*ooooi77i 
REAL  TITLECS)  00001774 

INTEGER  XC0DE«1RUn;nPAQE,LCT  00001779 

COMMON  /HEAD/  title, XC00E(VERS«lCvEL«0A7E,XRUN,HPA6CfNL06  0000177S 

COMMON/LINUM/LINE  00001777 

DATA  LCT/61/  0000177S 

C •..v.;...:.:«ooooi779 

LINEtLINE^N  000017S0 

XPdXNE.LT.LCT)  RETURN  000017S1 

LXNEvN«4  000017S2 

SO  HPAGEbnRACE^I  000017SS 

WRITE!*, 2030)  XCODE.XRUN, TITLE, VERS, LEVEL, date, mRAGE  000017S4 

2030  F0RMAT(«1*,IS,I*,9V,*AS,*  VERSION  *,29.1,*  f *, X*, *} *, 1 IX,  000017S9 

X A10,10X,*PA6E  *;IS/1X,'127(*»*))  000017S6 

RETURN  A 000017S7 

EnO  0000178S 


r>  o n o 


SUBROUTINE 

CRRX(N,NAME) 

00001789 

XBN  360 
VE9SXQN  C 

E.REXFENSTEIN 

LEVEL, 7204E1  . , 

fortran  XV  00001791 

0000179t 

XNTEOER  N 
REAL  NAME 

00001794 

00001799 

WRITE(6»t000)  NfNAMC  000017R7 

6000  F0RMAT(*0EXECUTI0N  7ERMXNATEn  DuC  TO  ERROR  RQi  IN  *»Aa)  00001798 

STOP  00001799 

END  OOOOISOO 


r 


SUBROUTINE  ERRM(N»NAME)  OOOOISOI 

*«*«*«************;****;;«**;**;*«**;«;*****;****«****««**ftt***ft**ooooisoi 

VERSION  1,0  LEVEt  rSOREl  OOOOISOS 

*********«***«*******************«******«******t*«******«**A****;*oooouo« 

integer  N OOOOlSOf 

REAL. NAHE  OOOOISOS 

NRITE(S»S100)  N,NAME  OOOOISOS 

SlOO  R0RHAT(*0ERR0R  no.  It, It,*  IN  *,A8/)  OOOOlSOf 

RETURN  OOOOISK 

END  OOOOlSlt 


c»  o o o o 


c 

10 

5010 


20 

30 

32 

fc032 

40 


50 


SUBROUTINE  INE(XC), 

IBM  3B0  E.REIFENSTeXN 

VERSION  1 LEVEL  720S02 

READS  AND  PRINTS  COMMENTS  CARDS 

REAL  name 

INTEGER  IP0RM,IP(3),C0M(13) tSLANK 

DATA  IF/IH  ,1H0,1HI/,NAME/SHINE/«BLANK/1H  f 

OOOOUl 

fortran  IV  ooooist 
ooootst 
ooooist 

OOOOISI 

ooooist 

00001St( 

READ(IC|5010)  IFORHfCQMfJP 

ooooist 

FORMAT (14X, A l,5Xtl2AA,A2,A2) 

0000  istj 

DO  20  III, 3 

0000  IStl 

IFdFORM.EO.IFd))  CO  TO  (30, SO, 40), I 

ooooist 

CONTINUE 

ooooist 

CALL  CRRX(20fNAME) 

ooooist 

CALL  LINESd),RCTURNS(32) 

ooooist 

WRITE(S,S032)  IFd),COH 

ooooist 

F0RMAT(A1,T21,12A4,A2) 

OOOOlSf 

60  TO  50 

OOOOISS 

CALL  PAGE 

OOOOISI 

Ii2 

OOOOISS 

60  TO  30 

OOOOISI 

IF(JF,NE. BLANK)  60  TO  10 

OOOOIBI 

RETURN 

OOOOISS! 

END 

000018S1 

U <J  u 


1 


SUBROUTINE  DAY  00001038 

***«******«*««*****ft******;*;*«*«**«*«;*****;************;********ooooi8ss 

VERSION  1.0  LEVEL  T0O9E1  OOOOISAO 

*******ft*****************ft****«ft******;*****ft*****************;*;*ooooi8ai 

real  TXTLE(S)  00001S4E 

COMMON  /HEAD/  TXTLE*ICOOC(yERS»LEveL,'DAT»ZRUN«NRACE*MLOC  OOOO1S0S 

c ••••.i..».,....«000018«4 

DAT  ■ DATE(O)  00001645 

RETURN  00001646 

END  00001647 


: i 


